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ABSTRACT 

We explore the collisional decay of disk mass and infrared emission in debris disks. With models, we show that 
the rate of the decay varies throughout the evolution of the disks, increasing its rate up to a certain point, which 
is followed by a leveling off to a slower value. The total disk mass falls off oc t~ 035 at its fastest point (where 
t is time) for our reference model, while the dust mass and its proxy - the infrared excess emission - fades 
significantly faster (ex r 8 ). These later level off to a decay rate of M tot (t) oc r 08 and M dust (t) or L lv (t) ex r 6 . 
This is slower than the oc t~ x decay given for all three system parameters by traditional analytic models. 
We also compile an extensive catalog of Spitzer and Herschel 24, 70, and 100 /im observations. Assuming a 
log-normal distribution of initial disk masses, we generate model population decay curves for the fraction of 
debris disk harboring stars observed at 24 /im and also model the distribution of measured excesses at the far-IR 
wavelengths (70-100 /im) at certain age regimes. We show general agreement at 24 /im between the decay of 
our numerical collisional population synthesis model and observations up to a Gyr. We associate offsets above 
a Gyr to stochastic events in a few select systems. We cannot fit the decay in the far infrared convincingly 
with grain strength properties appropriate for silicates, but those of water ice give fits more consistent with the 
observations. The oldest disks have a higher incidence of large excesses than predicted by the model; again, a 
plausible explanation is very late phases of high dynamical activity around a small number of stars. 
Finally, we constrain the variables of our numerical model by comparing the evolutionary trends generated 
from the exploration of the full parameter space to observations. Amongst other results, we show that erosive 
collisions are dominant in setting the timescale of the evolution and that planetesimals on the order of 100 km 
in diameter are necessary in the cascades for our population synthesis models to reproduce the observations. 
Keywords: methods: numerical - circumstellar matter - planetary systems - infrared: stars 



1. INTRODUCTION 

Planetary debris disks provide the most accessible means 
to explore the outer zones of planetary systems over their 
entire age range - from 10 Gyr to examples just emerging 
from the formation of the star and its planets at 10 Myr. De- 
bris disks are circumstellar rings of dust, rocks, and plan- 
etesimals, which become visible in scattered light and in- 
frared emission because of their large surface areas of dust. 
Because this dust clears quickly, it must be constantly re- 
plenished through collisions amongst the larger bodies, ini- 
tiated by the dynamical perturbing forces of nearby planets 
(|Wyattll2008|) . Thus, the presence of a debris disk signals 
not only that the star has a large population of planetesi- 
mals, but that there is likely to be at least one larger body 
to stir this population. The overall structures of these sys- 
tems are indicative of the processes expected to influence 
the structures of the planetary systems. Th ey result from 
sublimation temperatures and ice lines (e.g., Morales et al. 
20111) and sculpting by unseen planets (e.g.. iLiou & Zookl 
1999; Ouillen & Thorndike 2002; Kuchner & Holman 2003; 
Mora n et al.l 120041: Ilvloro-Martin et al.l 120051: IChiang et all 
20091) . as well as from conditions at the formation of the plan- 
etary system. 

However, debris disks un dergo significant evolution 
dRieke et al.ll2QQ5h I Wvattl 120081) . Studies of other aspects of 
disk behavior, such as dependence on metallicity or on bi- 
narity of the stars, generally are based on stars with a large 
range of ages, and thus the evolution must be taken into ac- 
count to reach reliable conclusions about the effects of these 



other parameters. Analytic models of the collisional processes 
within disks have given us a rough unde rstanding of their evo- 
lution (IWvatt et al] 120071: IWvattl 120081) . yielding decays for 
the steady state (constant rate of decay). Multiple observa- 
tional programs have characterized the decay of debris disks 
(e.g. ISpangler et alJ 120011: iGreaves & Wyattll2003l: iLiu et all 
[2003;|Rie ke et al. 2005; jMQoret al. 2006'; Sie gler et al 1)2001 
iGasparet al.1 [2009: Ca rpenter etakl 120091: iMoor et al.l I20TH) 
and indicate general agreement with these models. However, 
these comparisons are limited by small sample statistics, un- 
certainties in the stellar ages, and the difficulties in making 
a quantitative comparison between the observed incidence of 
excesses and the model predictions. 

In fact, mor e complex numerica l models of sin- 
gle systems (E iebault et ail 120031: lLohne etal.1 120081: 
Theb ault & Au gereau 2007; Gaspar et al.ll2012ah have shown 
that the decay is better described as a quasi steady state, with 
rates varying over time rather than the simple de cay slope of 
1 typic ally found in traditional analytic models. lLohne et all 
( 2008) present the evolution of debris disks around late-type 
stars, using their cascade model ACE. The y yield a total mass 
decay slope of 0.3-0.4. The models of Kenyon & Bromley 
(2008) yield a fractional infrared luminosity (fd = La ust /L*) 
de cay slope be t ween 0.6 and 0.8. The latest work presented 
by IWvatt et all (1201 1|) indicates an acceleration in dust mass 
decay, with the systems initially losing dust mass following a 
decay slope of 0.34, which steepens to 2.8 when Poynting- 
Robertson drag becomes dominant. For the same reasons as 
with the analytic models, these predictions are inadequately 
tested against the observations. We summarize the decay 
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slopes determined by observations and models in Table [T] 

In this paper, we compute the evolution of debris disk signa- 
tures in the mid- and far-in frared, using o ur numerical colli- 
sional cascade code dGaspar et al] 120 12a. Paper I hereafter). 
We examine in detail the dependence of the results on the 
model input parameters. We then convert the results into pre- 
dictions for observations of the infrared excesses using a pop- 
ulation synthesis routine. We compare these predictions with 
the observations; most of the results at 24 jim (722 solar- and 
359 early-type stars) are taken from the literature, but in the 
far infrared we have assembled a sample of 434 systems with 
archival data from Spitzer/MIPS at 70 fim and Herschel/PACS 
at 70 and/or 100 /mi. We have taken great care in estimating 
the ages of these stars. We find plausible model parameters 
that are consistent with the observations. This agreement de- 
pends on previously untested aspects of the material in debris 
disks, such as the tensile strengths of the pa rticles. Our ba- 
sic result confirms that of IWyatt et al] (l2007|) that the overall 
pattern of disk evolution is consistent with evolution from a 
log-normal initial distribution of disk masses. It adds the rigor 
of a detailed numerical cascade model and reaches additional 
specific conclusions about the placement of the disks and the 
properties of their dust. 

Although our models generally fit the observed evolution 
well, there is an excess of debris disks at ages greater than 
1 Gyr, such as HD 69830 and BD +20 207. We attribute 
these systems to late-phase dynamical shakeups in a small 
number of planetary systems. In support of this hypothe- 
sis, a number of these systems have infrared excesses dom- 
inated b y very small dust grains (identified by stro ng spectral 
features; iBeichman efafl 120051: ISong et al.ll2005h . The dust 
around these stars is almost certainly transient and must be 
replenished at a very high rate. For example, HD 69830 has 
been found to have three N eptune-mass planets within one 
AU of the star (lLovis et al.ll2006|) : they are probably stirring 
its planetesimal system vigorously. 

The paper is organized as follows. In section [2 we present 
the decay behavior of our reference model in three separate 
parameter spaces. In section [51 we introduce a set of carefully 
vetted observations that we will use to verify our model and 
to constrain its parameters, while in section [4] we establish 
a population synthesis routine and verify our model with the 
observed decay trends. In section [3 we constrain the param- 
eters of our collisional cascade model using the observations, 
and in section [6l we summarize our findings. We provide an 
extensive analysis of the dependence of the predicted decay 
pattern on the model parameters in the Appendix. 

2. NUMERICAL MODELING OF SINGLE DISK DECAY 

We begin by probing the general behavior of disk decay, us- 
ing the reference model pre sented in the second paper of our 
series dGaspar et al.ll2012b|) (Paper II, hereafter). Models fit- 
ted to the full set of observations will be discussed in Sections 
01 and We refer the reader to Papers I and II for the details 
of the model variables. We define the dust mass as the mass 
of all particles smaller than 1 cm in radius within the debris 
ring. In the Appendix, we analyze the dependence of the de- 
cay of a single disk on system variables also using the models 
presented in Paper II, and show the effects that changes in the 
model variables have on the evolution speed of the collisional 
cascade and/or its scaling in time. 

Our reference model (Paper II), is of a 2.5 AU wide 
(AR/R = 0.1) debris disk situated at 25 AU radial distance 
around an AO spectral-type star with a total initial mass of 1 



MEarth- The largest body in the system has a radius of 1000 
km. The dust mass-distribution of the model, once it reaches 
a quasi steady state, is well approximated by a power-law 
with a slope of 1.88 (3.65 in size space). In the following 
subsections, we describe the evolution of the decay of this 
model. We analyze the decay of three parameters: the total 
mass within the system, the dust mass within the system, and 
- to verify its decay similarity to that of the dust mass - the 
fractional 24 /im infrared emission (fd(24) = ^disk(24)/F*(24)). 

2.1. The decay of the total disk mass 

The decay of total disk mass is not observable, as a signif- 
icant portion of it is concentrated in the largest body /bodies 
in the systems, which do not emit effectively. As we show 
later, the evolution of the total mass is not strongly coupled to 
the evolution of the observable parameters, which is a "dou- 
ble edged sword". Fitting the evolution of the observables 
can be performed with fewer constraints; however, we learn 
less about the actual decrease of the system mass when us- 
ing a model that is less strict on including realistic physics 
at the high mass end of the collisional cascade. Also, the 
long-term evolution of the dust will be affected by the evolu- 
tion of the largest masses in a system, meaning that long-term 
predictions by models not taking this evolution into account 
correctly may be inaccurate. On the other hand, comparison 
between different collisional models and their collisional pre- 
scriptions is enabled by this decoupling. 

We show the evolution of the total disk mass of our refer- 
ence model in the top left and the evolution of the decay slope 
of the total mass in the bottom left panels of Figure [Q The 
evolution is slow up to 100 Myr (until the larger bodies settle 
in the quasi steady state), after which there is a relatively rapid 
decay. It reaches its steepest and quickest evolution around 1 
Gyr, when r ~ 0.35, where r is the time exponent of the de- 
cay (oc t~ T ). The decay then slows down; settling at r « 0.08. 
Such time dep endent variation in evolution speeds has also 
been shown by lLohne et al] (120081) . 

In the Appendix, we show that variations to the total initial 
disk mass only scale the decay trend in time (linearly), but 
not its evolution speed, meaning that more massive disks will 
reach the same r « 0.08, but at earlier times. Since our ref- 
erence model is a low-mass disk, the majority of observable 
disks will reach this slow evolutionary state well under a few 
Gyrs (a disk a hundred times more dense than our reference 
model will settle to its slow decay at « 1 Gyr). This property 
is used in the population synthesis calculations in section H 

2.2. The decay of the dust mass 

Analytic models of debris disks assume that they are in 
steady state equilibrium. Under such assumptions the dust 
mass decay is proportional to the decay of the total system 
mass. In reality, since there is no mass input at the high 
mass end, the systems evolve in a quasi steady state. Since 
mass evolves downwards to smaller scales within the mass- 
distribution, the further we move away from the high-mass 
cutoff, the better a steady- state approximation for the colli- 
sional cascade becomes. This is the reason steady- state ap- 
proximations for the observed decays have been relatively 
successful, but not exact. 

Our model shows a more realistic behavior. We show the 
evolution of the dust mass in the top middle, and the evolu- 
tion of the decay slope of the dust mass in the bottom middle 
panels of Figure Q] Since the final particle mass (size) distri- 
bution slope is steeper than the initial one (Paper II), dust mass 
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Table 1 

The decay trends in the literature, with proportionality of variables to time given as oc t~ T 






Paper 


M disk (t) f d (t) or f d{ 2A)(t) or M dust (f ) Exc (%) 


Notes 




Observations of ensembles of debris disks 



Silverstone (2000) 

Spangler et al. (2001) 

Greaves & Wvatt (2003) 

Liu et al. (2004) 

Rieke et al. (2005) 

Gasparetal. (2009) 

Moor et al. (2011) 

Soangler et al. (2001) 

Dominik & Decin (2003) 

Dominik & Decin (2003) 

Wvatt etal. (2007) 

Thebault et al. (2003) 

Lohne et al. (2008) 

Kenvon & Bromlev (2008) .... 

Wvatt etal. (2011) 

Wvatt etal. (2011) 

Wvatt etal. (2011) 

Wvatt etal. (2011) 

This work (valid for all systems) 
This work (valid for all systems) 



r = 1.75 
r = 1.76 

r <0.5* 

r = 0.7* 

T=1.0 

r = 0.3... 1.0 



r = 0.43 



Average f d fitted (clusters) 

Average f d fitted (clusters) 

Calculated from excess fractions assuming 

a constant distribution of dust masses 

Upper envelope of submm disk mass decay 

Spitzer MIPS [24] fraction 

Fitted published data between 10- 1000 Myr 

Dispersion between these extremes 



Analytic models of single debris disk evolution 



= 2.0 



T — 1.0 



= 2.0* 
= 2.0 
= 1.0 
= 1.0* 



Assumed steady-state 
Collision dominated removal 
PRD dominated removal 
Assumed steady-state 



Numerical models of single debris disk evolution 



r = 0.05 
r = 0.2 

r = 0.94 



r = 0.33 
r = 0.08 



r = 0.38* 
r = 0.3. ..0.4 
r = 0.6. ..0.8 

r = 0.34* 
r = 0.97* 
r = 2.8* 
r = 0.8* 
r = 0.6* 



Fitted between 3 and 10 Myr 



Above 100 Myr 
Below 200 Myr 
Above 2 Gyr 

PRD dominated above 10 Gyr 
At their fastest point in evolution 
At very late ages (quasi steady state) 



Population synthesis numerical models of debris disk evolution ^ 



This work (early types at 24 jum) 
This work (early types at 24 pm) 
This work (late types at 24 pm) 
This work (late types at 24 jum) 
This work (late types at 24 pm) 



T = 0.1 

r = 2.5 
r = 0.1 
r = 2.6 
r = 1.4 



10 - 250 Myr 
0.4 - 1 Gyr 
10-100 Myr 
0.2 - 0.4 Gyr 
0.6 - 10 Gyr 



* Decay timescale calculated for dust mass. 

t Disks placed at radial distances with disk mass distributions as described in Section [4] The decay describes the evolution of a disk population and not that 
of a single disk. 



will increase in the beginning of the evolution. The evolution 
speed increases up to around 0.01 Gyr, after which it stays 
roughly constant up to 0.1 Gyr. This is the period where the 
larger disk members settle into their respective quasi steady 
state. The evolution once again increases from 0.1 Gyr to a 
few Gyr, following the formation of the "bump" in the size 
distribution at larger sizes. The decay slows down again once 
the entire mass range has settled in its quasi steady state, with 
a decay ex t~° 6 . 

2.3. The decay of the fractional infrared emission 

Although our primary interest is the underlying mass and 
the largest planetesimals in a debris disk, the observable vari- 
able is the infrared emission of the smallest particles. The 
emission from the debris disk is calculated f ollowing algo- 
rithms similar to those in Gasparetal] (120081) . We show the 
evolution of the fractional 24 fim emission of our reference 
model in the top right, and the evolution of its decay slope in 
the bottom right panels of Figure [T] We follow the evolution 
of the fractional 24 fim emission instead of the fractional in- 
frared luminosity, as they will be identical in a quasi steady 
state and we spare integrating the total emission of the disk 
at each point in time. The plots clearly show that the evolu- 
tion of the emission is a proxy for the evolution of the dust 
mass in a system, as its decay properties mirror that of the 



dust mass. From hereon, we will only focus on the evolution 
of the infrared emission - which is the observable quantity - 
and neglect the dust mass. 

3. OBSERVATIONS 

We compiled an extensive catalog of 24-100 fim observa- 
tions of sources with reliable photometry and ages from var- 
ious sources. Spitzer 24 and 70 fim data fo r field stars wer e 
obtained from Sierchio et al. (in prep) and ISu et"aT] ((2006). 
We added 24 /im data from a number of stellar cluster studies 
(see Table El). Publicly avai lable PACS 70 and 100 urn data 
from the Herschel DEBRI S (Matthews 2008; Matthews^etaD 
[20l0h and DUNES (lEiroall2010l:lEiroa et alJl2010Ll201 lh sur- 
veys were also obtained from the HSA data archive. MIPS 24 
and 70 fim data for the stars in these surveys were also added 
to our analysis. 

The Herschel/PACS data were reduced using the Herschel 
Interactiv e Processing Environment (HIPE, V9.0 user release, 
lOttl 12010). We generated the calibrated Level 1 product by 
applying the standard processing steps (flagging of bad pix- 
els, flagging of saturated pixels, conversion of digital units to 
Volts, adding of pointing and time information, response cali- 
bration, flat fielding) and performed second-level deglitching 
with the "timeordered" option and a 20 sigma threshold to re- 
move glitches. This technique uses sigma-clipping of the out- 
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Figure 1. The decay of the reference model introduced in Paper II. The top row shows the decay of the total mass, the dust mass, and the fractional 24 fim 
infrared emission (fd(2A) = ^disk(24)/F*(24)); the bottom row shows the corresponding decay slopes for each parameter at the same points in time. The plots 
highlight that the decay is not a steady state process. 



lying flux values on each map pixel and is very effective for 
data with high coverage. After this stage the science frames 
were selected from the timeline by applying spacecraft- speed 
selection criteria (18'7s < speed < 20 // /s). The 1/f noise was 
removed using high-pass filtering with a filter width of 20 for 
the 100 fim data. This method is based on highpass median- 
window subtraction; thus the images might suffer from loss 
of flux after applying the filter. To avoid this we used a mask 
with 20 // radius at the position of our sources. After high-pass 
filtering we combined the frames belonging to the two dif- 
ferent scan directions and generated the final Level 2 maps 
using photproject also in HIPE. Aperture photometry was per- 
formed on the sources using a ^"radius, while the sky back- 
ground was determined with an aperture between 20" and 30". 
Six sub- sky apertures were placed within the nominal sky 
aperture with radii of 12", to estimate the variations in the 
sky background. Our self-calib ratio n of the data to determine 
the photospheric level (Section [33]) circumvents any residual 
calibration offsets. A summary of the photometry from the 
DUNES and the DEBRIS surveys as well as ages is presented 
in Table 12 The photometry of Sierchio et al. (in prep) will be 
published in an upcoming paper. 

3.1. Determining ages for the field sample stars 

Ages were estimated for these stars using a variety of in- 
dicators. Chromospheric activity, X-ray luminosity, and gy- 
rochronology as measures of stellar age are discussed by 
Mamajek & Hillenbrand (2008); we used their calibrations. 
To confirm the age estimates past 4 Gyr, we used a metallicity- 
corrected Mk vs. V - K HR diagram and found excellent cor- 
respondence between the assigned ages through extrapolation 
of their relations and the isochrone age. This work is dis- 
cussed in detail in J. Sierchio et al. (in prep). We also used val- 
ues of vsin/ > 10 km s -1 as indicators of youth, and logg < 4 
as an indicator of post-main-sequence status (when other in- 
dications of youth were absent). Our assigned ages and the 
sources of data that support them are listed in Table [2] We 
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Figure 2. The smoothed excess fraction decay curves at 24 fim for early- 
and late-type stars, with 1 a error bars. The late-type stars show a slightly 
quicker decay than the early-types. 

were not able to develop a rigid hierarchy among the methods 
in assigning ages, since occasionally an otherwise reliable in- 
dicator gives an answer that is clearly not reasonable for a 
given star - e.g., a low level of chromospheric activity can be 
indicated for a star whose position on the HR diagram is only 
compatible with a young age; HD 33564 is an example. 

3.2. The decay of planetary debris disk excesses at 24 jim 

Spitzer 24 /im data are reliable indicators of warm de- 
bris disk emission and have been used in a number of stud- 
ies (e.g.,[R ieke et alll2005l: iSu et al.ll2006l: ISiegler et al.ll2007l: 
Trilling et al. 2008; Gaspar et al. 2009). Given the uncertain- 
ties in the ages of field stars, stellar cluster studies, where nu- 
merous coeval systems can be observed, are strongly favored 
in disk evolution studies. The clusters included in our cur- 
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Table 3 

The excess fraction in [24] for early-type stars (A0-A9) and late-type stars (F5-G9) in clusters/associations. 



Name 


Age 




A0-A9 




F5-G9 


Excess Age 




[Myr] 


[#] 


[%] 


[#] 


[%] 


Reference 



P Pic MG 

LCC/UCL/US 

NGC 2547 

Tuc-Hor 

IC2391 

NGC2451B 

NGC2451A 

a Per 

Pleiades 

Hyades/Praesepe/Coma Ber 



123 

10-20 

30±5 

30±5 

50±5 

50±5 

65±15 

85* 

115±10 

600-800 



4/7 

42/89 

8/18 

2/5 

3/8 

0/3 

1/5 

5/26 
5/46 



c ? 1+14.9 
J/,i -18,0 

47.2+|-f 
444+11.7 

40.0+??! 

37.5! 

o +3b y 
+4.2 . 

20.0 + f 5 9 4 



1 9 2+9.9 
-§•3 
i o Q+ 6 - 3 
lu,y -3.0 



3/6 

42/92 

8/20 

0/1 

3/10 

6/16 

5/15 

2/13 

24/71 

1/47 



50.0±17.7 
45 7+5.3 

40.9±10.5 

+60 

30 +1 ^- 8 
DKJ - -1Q.Q 

37.5+^ 
333+f^ 
15.4!f 3 7 

33 - 8 ^50 

2 1 +46 



1 

3 

7 

1 

9 

11 

11 

13 

17 

20 



2 

4,5,6 

7 

8 

10 
12 
12 

14,15,16 
15,18,19 
21,22 



References. — (1) Rebull et al. (2008); (2) Ortega et al. (20®; (3) IChen et all l20TTh : (4) [Preibisch et al. (2002); (5) Fuchs et al. ( 20®; (6) 
Mamaiek et al. ( 2002); (7) Gorlova et al. ( 2007); (8) Rebull et al. ( 2008), with arbitrary errors adopted from similar age clusters; (9) Siegler et al. ( 2007); (10) 
Barrado v Navascues et al. (2004); (11) Balog et al. (2009); (12) Hunsch et al. (2003); (13) Carpenter et al. (2009); (14) Song et al. ( 2001); (15) Martin et al. 
( 2001) (16) Mamaiek & Hillenbrand ( 2008); (17) Sierchio et al. (2010); (18) Mevnet et al. (1993); (19) Stauffer et al. (1998); (20) Urban et al. (2012); (21) 
Gaspar et al. ( 2009); (22) Perryman et al. (1998) 



rent research (Table O have well defined ages and, more im- 
portantly, homogenic and reliable photometry. Unfortunately, 
getting an even coverage of ages using only clusters is not 
possible, especially for ages above a Gyr, which is why we 
combined the stellar cluster studies with field star samples. 
We include th e study of 24 /im excesses around early-type 
field stars by lSu et all (l2006|) . while the late-type stars are in- 
cluded from Sierchio et al. (in prep). We refer the reader to 
these papers on the methods used to determine precise ages 
and accurate photometry. We also include the Spitzer 24 /im 
measurements of the sources found in the DUNES and DE- 
BRIS Herschel surveys. Our final combined samples have 
722 and 359 sources in the late-type (F5-K9) and early-type 
(A0-F5) groups, respectively. 

At 24 /im, we determined excesses in the MIPS data by ap- 
plyin g an empirical rela tion between V - K and K- [24] (see, 
e.g., lUrbanet al.ll2012l) . We used 2MASS data for the near 
infrared magnitudes for many stars, but where these data are 
saturated we t ransformed he ritage photometry to the 2MASS 
system (e.g.. ICarpe nter 2001). In one case, we derived a K 
magnitude from COBE data, and in another we were forced 
to use the standard V - K color for the star, given its spectral 
type and B - V color (both of these cases are identified in Ta- 
ble [2]). We also determined an independent set of estimates of 
22 /im excesses from the WISE W3-W4 color. We found that 
on average this color is slightly offset from zero for stars of 
the spectral types in our study, so we applied a uniform cor- 
rection of -0.03. It is also important that the MIPS 24 /im and 
WISE W4 spectral bands are very similar, with a cuton filter 
at 20 and 19 microns, respectively, and the cutoff determined 
by the detector response (and with identical detector types). 
Not surprisingly, then, we found the two estimates of 24 /im 
excess to be very similar in most cases; where there were dis- 
crepancies, we investigated the photometry and rejected bad 
measurements. We then averaged the two determinations for 
all stars with measurements in both sets. We quote these av- 
erages, or the result of a single measurement if that is all that 
is available, in Table [2] 

For our current study, we are interested in the fraction of 
sources with excess as a function of stellar age. We defined 
a significant excess to occur when the excess ratio (defined 
as the ratio of the measured flux density to the flux density 
expected from the stellar photosphere) was > 1.1. Classi- 



cally, sources are binned into age bins and then the fraction 
of sources with excess is determined for each age bin. In- 
stead, we ran a Gaussian smoothing function over the ob- 
served age range, with a Gaussian smoothing width of 0.2 
dex in log(age). With this method, we generate smooth ex- 
cess fraction (defined as the fraction of the sample of stars 
with excess ratios above some threshold, in this case above 
1.1) decay curves. Errors of these decay curves were calcu- 
lated using the method described in Gasp ar et"aTI (|2009|) . Our 
final smoothed decay curves at 24 /im with ± 1 a errors for the 
early- and late-type stars are shown in Figure [2] The late-type 
stars show a slightly quicker decay between 0.1 and 1 Gyr, 
outside of the 1 a errors. We compare these decay curves 
with population synthesis models in Section HI 

3.3. The decay of planetary debris disk excesses at 
70-100 [im 

We measured excesses at 70 /im (MIPS) and 100 /im 
(PACS) relative to measurements at 24 /im (MIPS) and 12 /im 
(WISE W3), respectively. For the MIPS data, we computed 
the distribution of the ratio of 24 to 70 /im flux density, in 
units of the standard deviation of the 70 /im measurement 
(we rejected stars with 24 /im excesses in this distribution). 
The result, in the left panel of Figure [3 shows a well de- 
fined peak at the photospheric ratio. For the PACS data, we 
started by assessing the W3 data. The validity of these mea- 
surements was determined by a V-K vs. K- W3 color-color 
diagram, from which it was concluded that they are not signif- 
icantly degraded in accuracy due to saturation effects except 
for sources brighter than W3 = 2.7. For such sources, we de- 
termined an equivalent W3 magnitude from K measurements 
and V-K (i.e., from the same V-K vs. K- W3 diagram); they 
are noted in a footnote to the data table. We then computed 
the distribution of flux ratios (W3 vs. PACS 100 /im), also 
expressed in units of the quoted error at the far infrared wave- 
lengths. This distribution also has a prominent peak at the 
ratio corresponding to pure photospheric emission (see right 
panel of Figure 0. We fitted these two peaks with Gaussians 
between -4 and +2 standard deviations (we did not optimize 
the fit using larger positive deviations to avoid having it being 
influenced by stars with excesses). This procedure automati- 
cally calibrates the photospheric behavior, correcting for any 
overall departure from models, correcting any offsets in cal- 
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Figure 3. Determining the calibration between MIPS 24 and 70 /im, and WISE W3 and PACS 100 /im. For displaying the data, the bins were smoothed using a 
three-bin running average. 



ibration, and compensating for bandpass effects in the pho- 
tometry. We used these values to estimate the photospheric 
fluxes at 70 and 100 /im. In the case of 70 /im, we corrected 
the values for excesses at 24 /im by multiplying by the ex- 
cess ratio at this wavelength in all cases where it was 1.10 or 
larger. Smaller values are consistent with random errors and 
no correction was applied. The final result is two indepen- 
dent determinations of far infrared excesses. One is based on 
MIPS 24 and 70 /im measurements plus V-K and K. The sec- 
ond is based on WISE W3 measurements (with the exceptions 
noted in the footnote to the table) and PACS 100 /im ones. To 
test these results, we al so fitted stellar photospheric models 
(|Castelli & Ku rucz 2003) to the full set of photometry avail- 
able for each star from U through MIPS 24 /im and inspected 
the behavior of the MIPS and PACS 70 /im and PACS 100 /im 
relative to the photospheric levels predicted by these fits. This 
check neither called into question any of the excesses found 
previously, nor did it suggest additional stars with excesses. 

Based on the distributions in Figure [3 it appears that the 
quoted errors slightly underestimate the true ones, particularly 
for MIPS. We set the threshold for a detection to be 4 in units 
of the quoted measurement errors for MIPS and 3.5 for PACS. 
There are then 66 stars with far infrared excesses in our com- 
bined DEBRIS/DUNES sample, 57 of which have been mea- 
sured by both instruments. Of these 57, 37 have excesses at 
consistent and significant levels from both measurement sets. 
Eight stars are above the detection level for PACS but not for 
MIPS; five of these show positive signals in MIPS but below 
the threshold, and the remaining three are consistent with the 
PACS indicated excesses within two sigma. Twelve stars are 
above the detection threshold for MIPS, of which eleven show 
positive signals in PACS. That is, the measurements with the 
two instruments are very closely matched in their ability to 
identify far infrared excesses, and the measurements of the 
stars in common are very consistent. This is an important 
conclusion because not only are the measurements themselves 
independent, but our analysis of excesses maintained virtu- 
ally complete independence (with only the exception of PACS 
sources where we needed to supplement W3 with K data). 

By comparing the results from both instruments and main- 
taining their independence, we have been able to identify re- 
liably a set of stars with far infrared excess emission. We list 
the final photometric data for these sources in Table [2J For 



our current study we also include the 70 /im measurements of 
Sierchio et al. (in prep). Our final catalog of far-IR measure- 
ments of late-type stars totals 434 sources. 

4. POPULATION SYNTHESIS AND COMPARISON TO 
OBSERVATIONS 

In this section, we compare the decay of infrared excesses 
predicted by our model, using population synthesis, to the ob- 
served fraction of sources with excesses at 24 /im and to the 
distribution of excesses at 70- 100 /im. The two wavelength 
regimes are dealt with differently due to reasons explained in 
Section [3 

4.1. Disk locations 

By fitting black body e mission curves to IR S spectral en- 
ergy distributions (SEDs). lMorales et alJ (120111) found that the 
majority of debris disks have just a cold component or sepa- 
rate cold and warm components. Mostly independent of stel- 
lar spectral type, the respective temperatures for the warm and 
cold components yielded similar values. 

The warm component was found slightly above the ice 
evaporation temperature, with a characteristic temperature of 
190 K. While the systems around late-type stars have a nar- 
rower distribution in dust temperatures (99 to 200 K), the ones 
around A-type stars have a significantly wider one (98 to 324 
K). Assuming astronomical silicates as grain types in warm 
debris disks (where volatile elements are likely missing), we 
calculate the equilibrium temperatures of grains as a function 
of their sizes and radial distances around late- and early-type 
stars. We show these temperature curves in the top panels 
of Figure HI With green bands, we plot the particle size do- 
main that is most effective at emitting at 24 /im, when con- 
sidering a realistic particle size distribution within the disks 
( Gas par et al.ll2012b|) . With gray bands, we show the relative 
number of systems found by Morale set al.l (1201 ll) at various 
system temperatures. According to these plots, the most com- 
mon radial distance for warm debris disks (where the green 
band and gray bands intersect) is at « 4 AU around late-type 
stars. This can be seen in the figure because the temperature 
curves for 3.5 and 4.5 AU pass through the center of the in- 
tersection of the green and gray bands. A similar argument 
indicates a radial distance of ~ 1 1 AU for the early-type stars. 
However, a range of distances can be accommodated, espe- 
cially if one considers grains with varying optical properties. 
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Figure 4. Grain temperatures as a function of particle size, composition, radial distance, and the spectral type of the central star. The colored vertical bands yield 
the optimum particle size for emission around the certain systems, while the horizontal gray bands yield the relative number of systems found at each temperature 
(darker stand for more sources). The plots yield the general radial distance of warm and cold debris disks around different spectral type systems where the colored 
bands, the gray bands, and the temperature curves intersect. 



We performed similar analysis for the cold components, but 
only for the late-type sample, as we do not have a statisti- 
cally significant sample at old ages for the early-types. For 
the cold component analysis, however, we include a second 
grain-type, one that includes volatiles, as these disks are lo- 
cated outside of the snowline. We use the optical properties 
calculated by iMin et al.1 (1201 11) for a Si/FeS/C/Ice mixture, 
which have been used to successfully model the far-IR emis- 
sion a nd resolved imag es of Fomalhaut obtained with Her- 
schel (lAcke et al.l 12012). We show these plots (green band - 
astronomical silicates; red band - volatile mixture), in the bot- 
tom panels of Figure H The plots estimate the cold disks to 
be located at around 20 AU for an astronomical silicate com- 
position and around 25-30 AU for the volatile mixture. The 
latter estimate is more in agreement with the location of the 
Kuiper belt within our solar system. We can compare with 
disks around other stars by scaling their radii according to the 
thermal equilibrium distances, i.e., as (L*) 1 / 2 . The locations 
for grains of the ice mixture generally agree with these scaled 
radii. 



4.2. Modeling the 24 fim excess decay 

Based on the previous section, to model the decay of the 
warm components, we calculated the evolution of debris disks 
at radial distances between 2.5 and 10 AU with 0.5 AU incre- 
ments for late-type stars (GO), and at radial distances between 
9.0 and 14 AU with 1.0 AU increments for early-type stars 
(AO). The disk widths and heights were set to 10% of the 
disk radius, while the total disk mass was set to 100 M e . All 
other parameters were the same as for our reference model 
(Paper II). In Figure [5l we show the evolution of the model 
debris disk at 4.5 AU around a late-type star. The top left 
panel shows the evolution of the particle mass distribution in 
"mass/bin"-like units. The top right panel shows the evolu- 
tion of the SED of the debris disk, with the color/line cod- 
ing being the same as for the mass distributions. The SEDs 
were c alculated ass uming astronomical silicate optical prop- 
erties dDraine & Led 1 19841) . Both the mass distribution and 
the SED decay steadily in the even log- spaced time intervals 
we picked. The bottom left panel shows the evolution of the 
fractional 24 jim infrared emission, which (as with our ref- 



8 



GAspAr et al. 



i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — r 

yr 

100 yr 
■4000 yr 

0.16 Myr 

6.3 Myr 
"255 Myr 




Q II i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i 

-15 -10 -5 5 10 15 20 
log 10 (m/kg) 

i i| — i n-| — i 1 — r| — i — n — q — i 1 i| :i 1 — r| — i r-q — i nj- 



3 r 




1 r 




2 3 4 5 6 
log 10 [Age/yr] 



1 |a,m 1 (xm 1 00 (xm 

X 



1.0 
0.8 
0.6 
0.4 
0.2 h 



0.0 



■ 'i ■ ■ 'i 1 1 'i 




1 2 3 4 5 6 
log 10 [Age/yr] 



Figure 5. Evolution of the warm disk component model around a late-type star at 4.5 AU. Top-left panel: The evolution of the particle-mass distribution in 
"mass/bin"-like units. Top-right panel: Evolution of the SED of the disk (color coding is the same as for the top-left panel). Bottom-left panel: Evolution of the 
fractional 24 fim infrared emission as a function of age (the constant 5 Myr offset was applied later - see text). Bottom-right panel: The speed of the evolution of 
the fractional infrared emission. The evolution reaches its quickest point at 0.1 Myr, and settles to a oc t~° 6 evolution at around 10 Myr. The average mass disk 
in the population, which is barely detectable for a short period of time, would reach this at around 1000 Gyr; while a disk detectable between 0.1-1 Gyr reaches 
this quasi steady state around 1 Gyr. 



erence model in Section 3) shows varying speed in evolution. 
The color/line codes show the points in time that are displayed 
in the top panels. The speed of evolution is shown in the bot- 
tom right panel. The evolution speed curve is very similar to 
that of the reference model in Section |2 however, the evo- 
lution is much quicker. While our reference model settles to 
the oc t~°- 6 decay at around 100 Gyr, our warm disk model at 
4.5 AU already reaches this state at 10 Myr. There are two 
reasons for this behavior: 1) The disk evolves quicker closer 
to the star (the reference model was at 25 AU), and 2) the ex- 
tremely large initial disk mass (which was set to ensure cover- 
age at large disk masses as well) significantly accelerates the 
evolution. 

To compare these models with observations, we will use the 
excess fraction (fraction of a population with excesses above 
a threshold) as the metric, since this is the parameter most 
readily determined observationally. We calculate the fraction 
of sources with excesses at a given age using the decay of 
a single source and using a population synthesis routine, by 
making two assumptions: 

1. The distribution of initial disk masses follows a log- 
normal function. 



2. All systems initiate their collisional cascade at the same 
point in time during their evolution. This point can not 
be earlier than the time of planet formation. We fix t(0) 
at 5 Myr for our calculations. 

Both assumptions are plausible. Our first assumption is jus- 
tified by observati ons of protoplanetary disks in star form- 
ing regions. The Andrews & Williams (2005) study of the 
Taurus-Auriga star forming region showed that the protoplan- 
etary disk dust masses are distributed following a log-normal 
probability density function as 

r n 1 P / [ln(m)-/i] 2 | 
n(m\n,(i e )= -==Exp<^ — > , (1) 

where n(m) is the probability density of systems with initial 
masses of m, the "location parameter" of the log-normal dis- 
tribution is fx, and a e is the "scale parameter". We set the scale 
parameter to be equal to the widt h of the distribution of proto - 
planetary disk masses found by Andrews & Williams ( 2005), 
<Tg = 6.95 ± 0.06 (in natural log base). Since the peak in the 
mass distribution depends on the largest mass within the sys- 
tems and can be arbitrarily varied to large extent, the location 
parameter is found by fitting. We set the median (geometric 
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Figure 6. The best fitting mass population and its evolution for the warm component of late-type stars placed at 4.5 AU. Top panel: the fractional 24 /im emission 
decay curves (fd(24) = ^disk(24)/7 7 *(24)), shifted along the t~ l slope as a function of varying initial disk masses. The color code of the plot is proportional the 
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skewed. The black vertical line at fd(24) = 0.1 gives our detection threshold at 24 jim and the lower integration limit for our excess fraction decay calculations 



mean) of our log-normal distribution of masses to be equal to 

Cm = e /7 (2) 

where C is a scaling constant that yields the scaling offset 
between the median mass of the distribution and the mass of 
our reference model. 

The second assumption arises because the collisional cas- 
cades in debris disks cannot be maintained without larger 
planetary bodies shepherding and exciting the system. Ac- 
cording to core accretion models, giant pl anets such as Jupiter 
and Saturn form in less than 10 Myr (iPollack et al.l 11996k 
llda & Linl 120041) . while disk instabi lity models predict even 
shorter timescales (Boss 1997l l2Q0Tb . As planets form, simul - 
taneously, the protoplanetary disks fade (Hais ch et al.l l2001). 
and their remnants transition into cascading disk structures. 
Based on these arguments, our £(0) value of 5 Myr is rea- 
sonable. Our assumption ignores the possibility of later- 
generation debris disks. 

A useful property of collisional models is that their evolu- 
tion scales according to initial mass, which made the synthesis 
significantly simpler, as only a single model had to be calcu- 
lated. The flux / emitted by a model at time t with an initial 
mass m will be equal to a fiducial model's flux fo with initial 
mass mo at time to as 



/(0 = /oft>) — 
m 



(3) 



m 



(4) 



We verified that our model follows these scaling laws by run- 
ning multiple models with varying initial disk masses (see Ap- 
pendix). These relations are equivalent to a translation of the 
decay along a t~ x slope, which is why as long as the decay 
of single sources remains slower than t~ l , the decay curves 
will not cross each other. Similar behavior has been shown by 
lLohne et a l. ( 2008). This also means that each particular ob- 
served fit) value can be attributed to a particular initial disk 
mass and that at any given age the limiting mass can be cal- 
culated that yields a fractional infrared emission that is above 
our detection threshold. 

To compare with the observationally determined percent- 
age of sources above a given detection threshold, we need to 
find the initial mass whose theoretical decay curve yields an 
excess above this threshold as a function of system age. As 
detailed above, since the decay speed is always slower than 
t~ x , this will always be a single mass limit, without additional 
mass ranges. We can then calculate the cumulative distribu- 
tion function (CDF) of the log-normal function using these 
initial mass limits [m/(0] defined as 



CDF[m l (ty,ti,a e ] = - 



1+erf 




(5) 



10 



GAspAr et al. 




60 



lOMyr 



100 Myr 



1 Gyr 



lOGyr 



Age 



50 - 



40 - 



30 - 



- _ 20 - 





1 

R= 


: 10AU 




R= 


=11 AU 




R= 


12 AU ~~ 


- 


"hi,* 


*V - 












\ ♦ \ — 
\ * 
\ * 






\ * 






\ * ' 
\ * \ _ 












♦ 


Early type stars 


i 







10 Myr 



100 Myr 

Age 



1 Gyr 



Figure 7. The excess fraction decay curves calculated from our best fitting population synthesis models for warm disks at varying distances for the two different 
spectral groups. The right panel shows the models for the late-type stars, while the left panel shows them for the early-type stars. 



Although the distributions get skewed in the number den- 
sity vs. current mass (or fractional infrared emission) vs. age 
phase space, they remain log-normal in the number density 
vs. initial mass phase space, which is why this method can be 
used. The x 2 of our fitting procedure, where we only fit the 
location of the peak of the mass distribution, is then 



X 



^ {l-CDF[mi(t i y,^a e ]-F(t i )Y 



i 



(6) 



where Ffc) is the measured excess fraction at time tf, and 
dp (U) is the error of the measured excess fraction at time U. 
It is necessary to subtract the CDF from 1, because we are 
comparing the percentage of sources above our threshold and 
not below. 

In Figure [6l we show the best fitting mass population and 
its evolution for the warm component of late-type stars placed 
at 4.5 AU. The top panel shows the fractional infrared emis- 
sion decay curves, shifted along the t~ l slope as a function of 
varying initial disk masses. As the plot shows, the curves do 
not intersect, and they do not reach a common decay enve- 
lope (as is predicte d by analytic mod els that yield a uniform 
r 1 decay slope (e.g. Wyatt et al. 2007)). This also means that 
a maximum possible disk mass or fractional infrared emission 
at a given age, as predicted by these models, does not exist. 
The color code of the plot shows the number of systems at 
any given point in the phase space. While systems show a 
spread in fractional infrared emission up to ~ 100 Myr, after 
that they do increase in density along the decay curve of the 
average disk mass (shown with bold line) up to 10 Gyr, which 



is still faster (oc t 



-0.6 



...f 



-0.8 



) than the final quasi steady state 



decay speed of ex t~° 6 . The bottom panel shows the evolution 
of the number distribution as a function of fractional 24 fim 
emission at different ages (vertical cuts along the top panel). 
The initial distribution at age follows the initial mass dis- 
tribution's log-normal function; however, as the population 
evolves this gets significantly skewed. The black vertical line 
at fd(24) = 0.1 gives our detection threshold at 24 /im and the 
lower integration limit for our excess fraction decay calcula- 
tions. 

In Figure [71 we show the excess fraction decay curves cal- 



culated from our best fitting population synthesis models at 
varying distances for the two different spectral groups. The 
left panel shows the models for the late-type stars, while the 
right panel shows them for the early-type stars. The late-types 
can be adequately fit with models at 4.5 and 5.5 AU, which 
matches reasonably w ell to the temperature peak observed by 
Mor ales et al.l (|2011[) . Similarly, we get adequate fits to the 
early-type population with models placed at 1 1 AU, which 
is also in agreement with the temperature peak observed by 
iMorales et al.l (1201 1|) and our radial distance constraint. 

Our population synthesis routine yields excess fraction de- 
cays that are in agreement with the observations. This is the 
first time that a numerical collisional cascade code has been 
used together with a population synthesis routine to show 
agreement between the modeled and the observed decay of 
infrared excess emission originating from debris disks. The 
average disk mass predicted by our population synthesis has 
a total of 0.23 Mmooii, with a largest body radius of 1000 km. 
This yields dust masses of Md us t(< 1 cm) = 2.3 x 10~ 5 Mm 00 ii 
(^dust(< 1 mm) = 7.3 x 10~ 6 Mm oii). Our predicted aver- 
age dust mass is in agreement with the range of dust masses 
(2.8 x 10" 7 to 5.2 x 10" 3 M Mo0 n) observed bv iPlavchan et al.1 
(2009) for debris disks around young low-mass stars, deter- 
mined from infrared luminosities. 

4.3. Modeling thefar-IR (70-100 jim) excess decay 

According to section 14.11 to model the decay of the cold 
disks, we calculated the evolution of a disk placed at 15, 20, 
25, 30, and 35 AU around a late-type star. At these dis- 
tances, volatiles are a large part of the composition, which will 
change not only the optical properties of the smallest grains 
(see section |4J]), but also the tensile strength of the material. 
To account for this , we used the tensile strength properties 
of water-ice from Benz & Asphaug (199 91) and the eros ive 
cratering properties of ice from Koschny & Grtinl (l2001allb|) . 
For comparison, we repeated the calculations with the ten- 
sile strengths of basalt, as in our reference model. The emis- 
sion of the modeled particle size distributions was calculated 
assuming astronomical silicates for the reg ular basalt tensile 
strength models, and the volatile mixture (|Min et al.ll2011|) 
mentioned in section 14.11 for the water-ice tensile strength 
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models. 

Understanding and modeling the decay observed at far-IR 
wavelengths is significantly more difficult than it is for its 
shorter, 24 fim wavelength, counterpart. This is due to the 
non-uniform detection limits at longer wavelengths, which 
are frequently significantly above the stellar photospheric val- 
ues. Here, we will use the method developed by Sierchio et 
al. (in prep) to study the evolution of the far-IR excess, but 
slightly modified to use our calculated evolved fractional in- 
frared emission distributions. This new method quantifies the 
decay, taking into account both detections and non-detections 
and also the non-uniform detection limits. 

We define the significance of an observed excess as 



F-P _ Rf-1 



(7) 



where F is the detected flux, P is the predicted photospheric 
emission of the central star, while a is the error of the pho- 
tometry. We define Rf = F/P as the excess ratio of the source, 
and (Jr as the photosphere normalized error. As discussed in 
Section 3.3, a source is defined to have a detected excess if 
X > 3.5 (PACS/Herschel) or x > 4 (MlPS/Spitzer). 

The majority of the sources had both Spitzer 70 /im and 
Herschel PACS 100 /im data. We merged these data to simu- 
late a single dummy 85 /im datapoint as 
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Since the excess ratios at 70 and 100 jim are similar, when 
measurement was only available at a single band, it was as- 
signed to be at 85 /im. 

We separate our observed sources into three age bins that 
cover the age range between and 10 Gyr, the first bin in- 
cluding stars up to 1 Gyr (median age 0.34 Gyr), the second 
including stars with ages between 1 and 4 Gyr (median age 
2.6 Gyr), and the third with stars between 4 and 10 Gyr (me- 
dian age 6.4 Gyr). These age bins were chosen to include 
equal numbers of sources. 

We synthesize disk populations at 85 /im the same way as 
we did when modeling the 24 /im excess decay, assuming a 
log-normal initial mass distribution, with the scale parameter 
fixed at a 2 = 6.95, and varying only the location parameter of 
the distribution. 

Finally, we compare the calculated distribution at 0.34 Gyr, 
2.6 Gyr, and at 6.4 Gyr, to the observed first, second, and 
third data bin, respectively. Since the detection thresholds are 
non-uniform, instead of doing a straight comparison between 
the distributions, we calculate the number of possible detec- 
tions from our modeled distributions and compare with the 
observed distribution of excess significances (x's). Assuming 
that the model distribution does show the underlying distribu- 
tion of fractional far-IR excesses, we integrate the distribution 
upward from the detection threshold for each star in the cor- 
responding data bin. The detection threshold is given as 



7? — 1 

6 = 1 + 3- = 1+3^ — = 1 + 3(7/?. 
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(10) 



Integrating the distribution from the respective detection 
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Figure 8. The observed and modeled distribution of excesses at 25 AU 
around a late-type star, using water-ice particle tensile strength and a volatile 
mix for grain optical properties. The best fitting model using the fiducial 
basalt tensile strength and astronomical silicates for grain properties yielded 
similar distributions, only at 15 AU radial distance. 



threshold of each source yields the probability of detecting 
an excess at the given threshold according to the model. Sum- 
ming up these probabilities then yields the total number of 
predicted excesses that would be detected. This can then be 
compared to the actual number of observed excesses. The 
model that yields the best agreement for all three data bins 
consistently is defined as the best fitting model. 

In Figured we show the observed and modeled distribution 
of excesses at 25 AU, assuming water-ice tensile strength and 
the ice-mixture optical properties (the best fitting solution) in 
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15 20 25 30 35 15 20 25 30 35 



R (AU) R (AU) 

Figure 9. The relative difference between the predicted and observed number of far-IR excess sources at late-type stars, as a function of model radial distance. 
Left panel: Basalt tensile strength and astronomical silicate grain optical properties. Right panel: Water-ice tensile strength and volatile mixture grain optical 
properties. 



Table 4 

T he number of sources predicted vs. the number of sources observed in each age bin at different disk rad ii, assuming either astronomical silicates 
Draine & Lee (1984) for optical properties and basalt for grain tensile strength Benz & Asphaus (1999) or a Silicate/FeS/C-dust/Ice mixture 
IMin et al] (2011) for optical properties and ice for grain tensile strength Benz & Asphaug (1999). 



I N P /N for Silicates [and (g, (Basalt)] N P /N for Si/FeS/C/Ice Mixture [and 0£(lce)] 

R(AU) 10 Myr < t < 1 Gyr 1 Gyr < t < 5 Gyr 5Gyr<t<10Gyr 10 Myr < t < 1 Gyr lGyr<t<5Gyr 5 Gyr < t < 10 Gyr 



15 


32.35/33 


26.07/26 


15.36/14 


37.65/33 


16.11/26 


9.92/14 


20 


31.09/33 


29.21/26 


24.80/14 


34.14/33 


22.72/26 


10.96/14 


25 


24.29/33 


25.47/26 


24.14/14 


32.36/33 


26.82/26 


15.55/14 


30 


28.02/33 


30.60/26 


29.85/14 


29.12/33 


27.15/26 


20.94/14 


35 


25.42/33 


28.87/26 


28.79/14 


26.94/33 


27.83/26 


25.19/14 



the three separate age bins. The observed sources are com- 
pleteness corrected (by using the error distribution of the data 
itself), and sources below Rf < 1 are not shown. The pan- 
els display the number of sources observed and predicted by 
our calculations in each given age bin. We emphasize, that 
the numbers of predicted sources are not determined based 
on these binned emission plots, but with the method detailed 
above. These plots show the emission distribution predicted 
by our fits and compares it with the completeness corrected 
observed distributions. The distributions are scaled to the total 
number of sources. The best fit for the basalt tensile strength 
and astronomical silicate optical property model (which looks 
almost identical to the ice mixture/strength solution plotted) 
was at ~ 15 AU, which is clearly inwards of the predictions 
we made in section I47T1 and inwards of the cold disk compo- 
nent of our solar system. However, the water-ice composition 
and tensile strength model yields a fit that is in agreement 
with the predictions. In Table 01 we tabulate the number of 
predicted and observed sources for both models at various ra- 
dial distances, and in Figured we plot the relative differences 
between these numbers and show the predicted radial location 
of the disks with a red band. 

4.4. Disk incidence for old stars 

At 24 /im, our model suggests there should be virtually no 
detected debris disks around stars older than 1 Gyr. Nonethe- 
less, there are a number of examples, and examination of their 
ages indicates that they are of high weight. This resu lt im- 
plies that the simple assumption (e.g.. lWyatt et al.ll2007l) that 
debris disks can be modeled consistently starting from a log- 



normal initial mass distribution is successful up to about a 
Gyr, but that there are additional systems around older stars 
above the predictions of the simple model. We attribute these 
systems to late-phase dynamical activity that has led to sub- 
stantial enhancements in dust production. Two examples are 
HD 69830 and BD 20 207. Both of these systems have strong 
features in their infrared spectra that indicate the emission 
is dominated by small grains that must b e recently produced 
(iBeichman et al.ll2005l: ISong et al.ll2005h . which supports the 
hypothesis that they are the sites of recent major collisional 
events. Similarly, although our model successfully matches 
the numbers of detected disks in the far infrared, the obser- 
vations find many more large excesses than predicted (Figure 
[51 bottom panel). A plausible explanation would be that the 
outer, cold disk component can also have renaissance of dust 
production due to late phase dynamical activity. 

5. CONSTRAINING MODEL PARAMETERS WITH OBSERVATIONS 

We ran more than a hundred extra models, taking our best fit 
to the decay of the warm component of late-type debris disks 
as the basis, to test the dependence of the decay on the vari- 
ables of the model. We varied each model parameter within 
a range of values and performed the same population synthe- 
sis routine and fitting as we did in Section 4. Of these, nine 
variables show signs of having some effect on the evolution of 
the excess fraction decay curve. In Figure [lOl we present the 
reduced x 2 minima at each value of these nine parameters. 

Variables a and b of the cratered mass equation had the 
strongest effect on the slope of the evolution (see Appendix) 
and also strongly affect the population synthesis fits. Values 
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Figure 10. The reduced x 2 minima of the model population fits for each tested value of the selected nine model variables that have the largest effects within the 
fits. Red lines show the values of variables used in the base model. 



of a and b that describe materials that are softer in erosive col- 
lisions (a > 10~ 5 J kg -1 , b > 1.27) can be generally ruled out 
by our analysis for the warm component of debris disks. Our 
analysis also shows that the measured values of these vari- 
ables, which we used in our reference models, yield accept- 
able fits with our population synthesis routine. This is similar 
to the effect we observed when using water-ice erosive prop- 
erties for the cold disk components in the previous section. 

While the value of the slope of the tensile strength curve s 
significantly affects the slope of the particle mass- distribution 
(lO'Brien & Greenb erg 2003: lGaspar et al.ll2012bl) . it does not 
affect the decay of the fractional infrared emission to the level 
where we would observe offsets between the modeled and ob- 



served rates. However, we do have a best fit at its nominal 
value. 

The effects of varying S and Q s are roughly the same as 
when varying s. As it turns out, the exact value of the ten- 
sile strength law does not strongly influence the decay of the 
excess fractions in a population of debris disks. However, 
choosing a higher value for fM, which gives the interpola- 
tion distance between the erosive and catastrophic collisional 
domains, does result in less acceptable fits. This is an ar- 
bitrarily chosen numerical constant, and this analysis shows 
that choosing its value wisely is important. Based on these 
findings, we conclude that for our cold disk models in the 
previous section, the changes to a and b when assuming a 
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water-ice strength for the erosive collisions had a larger effect 
on the evolution than the changes to the catastrophic collision 
properties of the tensile strength curve. 

While varying r] (the initial particle mass distribution slope) 
of a single disk will have significant effects on the timescale 
of its evolution (see Appendix), it does not strongly determine 
the timescale of the excess fraction evolution of a population. 
To compensate for the offset in timescales, the average disk 
mass varies from population to population (within an order 
of magnitude). Testing the actual value of the initial particle 
mass distribution is possible, by comparing the disk mass dis- 
tributions predicted for each population to observations (such 
as in young clusters). 

Varying the maximum mass m max of the system did not have 
a large effect on the population synthesis fits above 10 18 kg, 
which reinforces our previous statement that it is the dust 
density of the model that matters and not an absolute total 
mass or largest mass in the system, which are redundant vari- 
ables. However, very low maximum mass systems (< 10 18 kg 
- « 100 km diameter) will result in decays that are inconsis- 
tent with our observations. This also has the important con- 
sequence, that the evolution of the planetary systems has to 
reach the point where bodies on this size scale are common in 
order to have a "successful" collisional cascade. 

The radial distance of th e mo del (R) obviously is the dom- 
inant parameter. In section l4~2l we showed that the best fit of 
our model to the observations is at R w 4.5 - 5.5 AU, which 
agrees with the thermal location predicted by Morale s~et al.l 
( 201 lj). Here, we show the quality of the fits when varying 
the radial distance between 2.5 and 10 AU. Placing the disks 
closer than 4 or further than 8 AU yields a population de- 
cay that is inconsistent with the observations. This value can 
likely be modified to some extent by varying some of the other 
input variables of the model. 

6. CONCLUSIONS 

In this paper, we present a theoretical study of the evolution 
of debris disks, following their total disk mass (M tota i), dust 
mass (Mdust)» an d fractional 24 /im infrared emission (fd(2A))- 
We use the numerical code presented in Paper I that models 
the cascade of particle fragmentation in collision dominated 
debris disk rings. 

Observational studies in the past decades have shown that 
the occurrence and strength of debris disk signatures fade 
with stellar age (e.g ISpangler et al]l200U iRieke et alJl2005t 
iTrilling et al.l 120081: ICarpenter et al.ll20Q9h . Analytic models 
of these decays explained them as a result of a steady-state 
(equilibrated) c ollisional cascade between the fra gments (e.g., 
ISpangler et all 120011: iDominik & Decinl 120031: IWvatt et al.l 
2007), which results in a decay timescale proportional to ex t~ l 
for all model variables (M tota i, Md US t> fd(24))> Analysis of the 
observed decays of stellar populations, however, has shown 
that the dust mass and the fractional infrar ed emission - the 
observable parameters - decay less quick ly (iGreaves & W vatt 
12001 Liu etal.ll200l iMoor et al.l I20TTL e.g.,). Slower de- 
cays have also be en modeled by comp l ete numerical cas- 
cade models (e.g.. [Thebault et al.l 120031: lLohne et al.l 120081: 
Kenyon & Bromley 2008|). Numerical codes yield slower de- 
cays because they model the systems as relaxing in a quasi 
steady state, instead of in complete equilibrium. This means 
that mass is not entered at the high mass end into the sys- 
tem (like in an analytic model), but is rather conserved. The 
remaining discrepancies among the numerical models are re- 
sults of the different collisional physics and processes mod- 



eled within them. 1 

Our calculations show that the evolution speed constantly 
varies over time and cannot be described by a single value. 
Since the fractional infrared emission is a proxy for the dust 
mass, their decays closely follow each other. At its fastest 
point in evolution, the total mass of our models decays as 
Motai oc r 33 , while the dust mass and fractional infrared 
emission of the single disk decays oc t~° mS . At later stages 
in evolution these slow down to oc t~ om and ex t~° 6 , respec- 
ti yely. These results are mos tly in agreement with the models 
of lKenyon & Bromleyl(|2008|) . We roughly agree w ith the dust 
mass decay predicted by the IWvatt et aD (|2011[) models up 
to the point where Poynting-Robertson drag (PRD) becomes 
dominant in their models (although their models decay some- 
what faster than ours, possibly due to the constant effects of 
PRD). 

We perform a population synthesis routine, assuming a log- 
normal probability distribution of initial disk masses. We cal- 
culate excess fraction decay curves, which we fit to the ob- 
served fraction of warm debris disks at a 10% excess thresh- 
old at 24 /im. Our fits show a good agreement between the 
calculated and observed decay rate of the fraction of debris 
disk sources around both late and early-type stars, with mass 
ranges in agre ement with the distributio n of protoplanetary 
disk masses ([Andrews & Williams 2005). We also analyze 
data from the MlPS/Spitzer and the DEBRIS and DUNES 
Herschel Space Observatory surveys. Taking into account 
the non-uniform detection thresholds at these longer wave- 
lengths, we also show roughly good agreement between the 
number of sources predicted to have an excess from our pop- 
ulation synthesis routines and that observed within these sur- 
veys. 

There are a small number of bright debris disks at 24 /im 
around old stars that are not predicted by the simple decay 
from a log-normal starting distribution; they [HD 109085, 
HIP 7978 (HD 10647), HIP28103 (HD 40136), HIP40693 
(HD 69830)] probably represent late-phase dynamical activ- 
ity. Similarly, the model fails to fit the large excesses in the 
far infrared around old stars, again consistent with late-phase 
activity around a small number of stars. 
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log is funded by the Deutsches Zentrum fur Luf- und Raum- 
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APPENDIX 

A. THE SYSTEM VARIABLES AND THEIR EFFECTS ON THE EVOLUTION OF THE COLLISIONAL CASCADE 

As we have shown in Section 5, varying the parameters of the model can affect the results of the population synthesis. Here, 
we analyze the effects of varying them on a single system. We summarize and describe the variables of the model in Tabled 

A.l. Evolution of the system mass 

We show the total mass decay curves as a function of model variables in Figure [TT] and the evolution of the power exponent of 
time in the decay of the total mass [M(t) ex t~ T ] as a function of these collisional variables in Figure[l2l The figures include plots 
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Table 5 

Numerical, Collisional, and System parameters of our model and their fiducial values 



Variable Description Fiducial value 



System variables 

~p Bulk density of particles 2.7 g cm _j 

m m i n Mass of the smallest particles in the system 1.42x 10 -21 kg 

m max Mass of the largest particles in the system 1 . 13 x 10 22 kg 

Mtot Total mass within the debris ring 1 M§ 

rjo Initial power-law distribution of particle masses 1.87 

R Distance of the debris ring from the star 25 AU 

AR Width of the debris ring 2.5 AU 

h Height of the debris ring - 2.5 AU 

Sp Spectral-type of the star A0_ 

Collisional variables 

7 Redistribution power-law 11/6 

fix Power exponent in X particle equation 1.24 

a Scaling constant in M cr 2.7 x 10~ 6 

b Power-law exponent in M cr equation 1.23 

/m Interpolation boundary for erosive collisions 10 -4 

f Y Fraction of Y/M cr 0.2 

/™ ax Largest fraction of Y/X at super catastrophic collision boundary ... 0.5 

<2sc Total scaling of the Q* strength curve 1 

S Scaling of the strength regime of the Q* strength curve 3.5 x 10 7 erg/g 

G Scaling of the gravity regime of the Q* strength curve 0.3 erg cm 3 /g 2 

s Power exponent of the strength regime of the Q* strength curve . . . -0.38 

g Power exponent of the gravity regime of the Q* strength curve 1.36 

Numerical parameters 

S Neighboring grid point mass ratio - 1.104 

Constant in smoothing weight for large-mass collisional probability 10 6 m max 

P Exponent in smoothing weight for large-mass collisional probability 16 



for the twelve variables that have the largest effect on the evolution, out o f the total twenty-four variables (see Paper I). These 
decays are compared to that of our r eference model, detailed in section [27TT 
In our code, we use the models of Benz & Asphaug (1999) to estimate the collision tensile strengths of particles, written as 

Q * (a) = i -*JA^Q\ S (^Y +G J^Y] , (Ai) 

erg kg L VI cm/ VI cm/ J 

where a is the target particle's radius, Q sc is the total scaling of the curve, S is the scaling of the curve in the "strength dominated" 
regime, s is the power exponent of the target radii in the "strength dominated" regime, G is the scaling of the curve in the "gravity 
dominated" regime, p is the bulk density of the particles, and g is the power exponent of the target radii in the "gravity dominated" 
regime. Of these, we show the effects of varying Q sc , G, and g, as varying S and s will not have a significant effect on the decay of 
the total mass, because they influence the low mass end of the distribution. Increasing or decreasing the total scaling will speed 
up the evolution of the total mass. Decreasing the total scaling of the tensile strength curve will soften the materials, resulting in 
a faster decay. Increasing it, however, will strengthen the materials, which will make the largest bodies "indestructible", resulting 
in a faster decay in the number of bodies just below the high mass end. A similar effect can be seen when G is varied. 

The total mas s crater ed in an erosive collision is calculated in our model by applying the experimental results of 
Koschnv & Grim (2001a : b). This mass is a function of a (scaling constant) and the projectile's energy to the b power. Vari- 
ations in these constants will affect how quickly the largest bodies erode and subsequently, the evolution of the total disk mass. 
When softer material properties are used (a and b increase), the decay is quicker, for example meaning debris disks composed of 
ice are likely to disappear in a sh orter timescale than rocky debris disks. 

The I Koschnv & Griinl (l2001bl) formula for cratered mass in an erosive collision is only valid for relatively small cratered 
masses. The cratered mass given by the formula can exceed M/2 even below the Erosive/Catastrophic collision boundary. We 
thus interpolate the cratered mass from /m = M cr /M to the boundary vi a methods given in Paper I. Assigning it a very small value 
basically eliminates the erosive formula of iKoschny & Griinl (l2001bl) and uses an interpolative formula for the entire domain. 
However, a larger value is likely to overestimate the cratered mass in an erosive collision near the erosive/catastrophic collision 
boundary. Our approach was to use a conservative value within these extremes. 

The number densities of fragments created in collisions in our model follow a power-law distribution. The slope of this 
distribution is given by 7, and only very minor effects can be seen when varying its value. The actual redistribution function 
has been a long researched topic within collisional systems, with some research s howing that double or even triple power-law 
functions are the best to describe the fragment distributions dDavis & Ryan 1990). According to our models, as long as the 
distribution function is within reasonable limits (7 < 1.99 - mass is concentrated in the largest fragments), there is not much 
difference in the decay of the total mass when varying its value. 

The total mass within the disk, M toU sets the scaling of the particle size distribution (as do m max and the volume of the disk). 
When scaling the initial total mass in the system, with all other parameters fixed, the evolution of the total mass is shifted in time, 
with the systems reaching their points of fastest decay at later points in time. This property is used in our population synthesis 
calculations in section 01 
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Figure 11. Evolution of the total disk mass as a function of selected parameters that have the largest effect on the timescale of the evolution. The numerical 
variable p modifies the smoothing function of the collisional cross section of the largest bodies in the system. The smoothing function only varies the evolution 
of the total mass (shown here), but does not affect the evolution of the dust mass or the fractional infrared emission. 
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Figure 12. The evolution of the power exponent of time in the decay of the total system mass [M(t) oc t~ T ] as a function model variables for decay curves shown 
in Figure [TT] At its fastest point, our reference model decays with r = 0.33, while all models reach a fastest point between 0.3 < r < 0.4. The evolution of the 
power exponent is characteristic for all models, with an acceleration in evolution up to a certain point, from whereon the evolution of the total mass slows down. 
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The decay is dependent on the mass of the largest body m max , which is usually arbitrarily chosen in the numerical models. This 
shows in our calculations in Section [3 where going below a largest body mass of 10 18 kg 100 km diameter) will result in 
decays that are inconsistent with our observations. When testing this for a single system, we set the total mass of the system to a 
value that yielded the same scaling of particle densities as the fiducial model had. This way we guaranteed that our calculations 
were only testing how varying the cutoff of the mass distribution affects the evolution. 

The slope of the initial distribution, rj, determines the number of dust particles when the collisional cascade is initiated. Our 
convergence tests (Paper I) have shown that the systems will reach collisional equilibrium from all initial distribution slope values. 
However, the time when the system reaches equilibrium will depend on the value of rj. A system will be able to reach equilibrium 
from slope values lower than the steady- state cascade distribution faster than from steeper slopes, as it is easier to produce and 
build up dust sizes, than to remove the large massive particles from the highest masses. 

One of the most important system variables is R, the distance of the disk from the central star. This parameter has many effects, 
as it sets the collisional velocity, thus the collisional energy of the particles and their collisional rate. It will also set the removal 
timescale for the blowout particles and is a variable in the volume of the disk, thus it sets the number density of the particles in 
the disk as well. Increasing the the radial distance will decrease the evolution rate of the disks, as shown in our Figures [TT] and 
IT2l with the fastest evolution setting in at later points in time. 

The last parameter we analyze is p, which is a variable that sets the smoothing function of the collisional rates for the largest 
bodies in the system (Paper I). Its value only affects the evolution of the largest masses, thus also the evolution of the total mass 
in the system. 

A.2. Evolution of the dust mass and fractional infrared emission 

As we have shown before, the fractional infrared emission is a proxy of the dust mass in the system, meaning the decay curves 
and the analysis we give for the fractional infrared emission are generally identical to the one we would give for the dust mass in 
the system. For said reasons, we omit the plots for the evolution of the dust mass. 

The emission of the particles depends on their temperat ures, their sizes, and mater ial and wavelength dependent optical prop- 
erties, such as their absorption coefficients. We assume a Castelli & Kurucz (2003) intensity emission model for the stars and 
astronomical silicate optical constants for the particles ( Draine & Lee 1984), when calculating their equilibrium temperatures and 
emission. 

We analyze the same parameters as in the previous subsection, with the exception of G, g, and p, which are replaced by S, 
s, and S. In Figure [131 we show the decay of the fractional infrared emission as a function of the model variables that have the 
largest effect on it, while in Figure [141 we show the power exponent of time in the decay. These figures can be compared the 
the evolution of th e inf rared emissions of our reference model, which is plotted with a thick solid line in the Figures and also 
analyzed in section l23l 

The variables of the tensile strength curve that determine the strengths of the gravity dominated larger bodies (G, g) do not 
affect the evolution of the dust distribution, while the variables that determine the tensile strengths of the smaller particles (S, 
s) do. Increasing or decreasing the scaling of the tensile strength law (Q sc ) increases the evolution speed for the dust mass, and 
thus the fractional infrared emission. At increased material strengths the quick decay of the largest bodies affects the evolution 
of the dust mass, while for softer materials a general faster decay of the entire distribution can be seen (see Figure 4 in Paper 
II). However, only significant decreases in the strength scaling S will have noticeable effects in the evolution of the fractional 
infrared emission. Increasing the steepness of the tensile strength law s will shift the evolution in time. Of all collisional variables, 
arguably b and a are the most important. As expected, using softer erosive material properties (larger b and a) speeds up the 
evolution of the dust mass (and with that the fractional infrared emission). 

Changes in f M and 7 affect the evolution of the fractional infrared emission similarly to that of the total mass. Increasing the 
largest body in the system (m max ) slows down the evolution of the collisional cascade, with models reaching their peak dust mass 
evolution at later stages, while increasing the total mass m tot in the system will speed the evolution of the system, with higher 
total mass systems reaching their peak evolutionary point earlier on. Systems initiating their collisional cascades with varying 
initial mass-distribution slopes (7/0) will reach their quasi steady state dust mass decay (the peak of evolution speed) roughly at 
the same time, even though the beginning of the evolution is dependent on it. Debris rings located at different radial distances (R) 
will evolve with speeds associated with their orbital velocities, shifting the onset of their quasi steady state decay to later points 
in time for disks at larger radial distances. 

Since p is the smoothing function of the largest bodies, it also does not affect the evolution of the dust mass; however, the 
neighboring grid point mass ratio (5) will be numerically important. In Figure[l3l we show that our models converge in dust mass 
decay at around 400-800 grid points, while using a less dense grid will result in numerical errors. 

A. 3. Conclusion 

Our analysis above has revealed that erosive collisions are dominant in shaping the evolution of a debris disk. The evolution 
speed of our model is determined primarily by the variables (a and b) of the cratered mass equation, when considering fixed 
system variables. This is not that surprising, considering that b also was found to be dominant in determining the mass-distribution 
slope (Paper II), and that our population synthesis analysis in section \5\ also revealed that our fits are sensitive to the values of 
a and b. The evolution is much less dependent on the catastrophic t ensile strength than on the e rosive, which is surprising, 
considering the dependence of the particle mass-distribution slope on s ( O'Brien & Greenberg 2003, Paper II). 

The measurements of iKoschny & Gr iin (|2001alb|) give the value o f a fo r silicates as 2.7 x 10~ 6 kg J -1 , and 6.2 x 10~ 5 kg J -1 
for ice; and a val ue of 1.23 for b. Measure ments by lHiraoka et alJ (120081) yield a b value of 1.15, which is in agreement with 
the value given by Koschny & Griin ( 200 la b) and yields an even better fit for our population synthesis constraint in section[5](a 
values cannot be compared as the papers used slightly different equations). 
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Of the system variables, the evolution will most strongly depend on 770, R, and m max . The evolution converges above m max = 
1 x 10 19 kg 200 km diameter), which most systems likely achieve (considering the asteroid sizes in our Main Asteroid Belt), 
making this variable less important for realistic conditions. Although 770 is difficult to constrain, it is likely that the system will 
form with a mass-distribution slope with a value close to its quasi steady state solution. However, even if they do not, they still 
can adequately reproduce our observations according to our population synthesis calculations in section [3 

The radial distance of the disk is the overall dominant parameter in determining the evolution of a single disk, when all realistic 
conditions are considered. It influences the evolution by three independent effects, all acting in the same direction. At larger 
radii, the collisional velocity will be lower (thus the collisional energy will be lower), which lowers the effective mass range a 
particle can interact with. The reduced collisional velocity also reduces the collisional rate. Finally, an increase in radial distance 
increases the effective volume the disk encompasses (for the same amount of mass and disk aspect ratio), also resulting in reduced 
collisional rates. 
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Figure 13. Evolution of the fractional 24 /im infrared emission (fd(24) = ^disk(24)/F*(24)) a function of selected parameters that have the largest effect on the 
timescale of its evolution. The characteristic bump seen in the evolution of the dust mass is reflected in the evolution of the infrared emission as well. The bump is 
followed by a drop in emission, which follows the same power-law as the drop in dust mass. Systems generally reach the quasi steady state decay at ~ 100 Myr, 
although variations in this are seen as a function of model parameters. 
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Figure 1 4. T he evolution of the power exponent of time in the decay of the fractional infrared emission as a function model variables for decay curves shown 
in Figure[l3] At its fastest point, the infrared emission our reference model decays with r = 0.8, while all models reach a fastest point between 0.6 < r < 0.9. 
These model results generally agree with observations of disk decay. 
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1.1Z 


1 no 

l.UZ 


y. 10 


a aa 
0.00 


n /io 

-U.4Z 


-5.07 


4.68 


-2.36 


N 


3 


2;4;15;16 


HDO 16970 


A3V 


0.50 


368 


40.91 


19.81 


361 


3.62 


0.98 


39.50 


5.23 


-0.27 


18.20 


5.76 


-0.28 


N 


2 


7 


HDO 17051 


F8V 


1.20 


163 


19.35 


8.84 


171 


1. / 1 


1.05 


01 O/l 
zi.y4 


J.UD 


n q< 
U.5D 


15.27 


3.30 


1.95 


N 


3 


2;3;4;5;6;14 


HD0 17093 


A7III 


0.58 


100 


11.45 


5.56 






1 C\A 
1.U4 








14.43 


3.10 


2.86 


N 


2 


7 


HDO 17206 


F7V 


0.75 


342 


38.72 


19.16 


346 


3.46 


1.01 








9.61 


3.59 


-2.66 


N 


3 


2;3;4;16 


HDO 18978 


A3IV-V 


0.70 


234 


26.79 


13.07 






1.03 








22.90 


3.60 


2.73 


N 


2 


7 


HD019107 


A8V 


0.07 


87 


9.92 


4.82 






1.03 








15.62 


2.96 


3.65 


? 


2 


7 


HDO 19305 


K5 


2.55 


45 


4.97 


2.32 


44 


0.45 


0.97 


-16.16 


5.62 


-3.76 


3.10 


2.43 


0.32 


N 


1 


I 


HD020010"f 


F6V 


4.80 


663 


71.31 


44.25 


630 


6.30 


0.95 


102.00 


5.21 


5.89 


36.52 


2.73 


-2.83 


N 


2 


3;21 


HD020320 


A9mA9V 


0.80 


153 


16.65 


8.16 


160 


1.61 


1.05 








99.90^ 


4.02 


22.82 


Y 


2 


7 


HD020630 


G5Vv 


0.40 


352 


41.44 


18.05 


366 


3.66 


1.04 


35.64 


6.91 


-0.84 


25.16 


3.76 


1.89 


N 


3 


2;4;6;12;13;14 


HD020794 


G8V 


6.20 


722 


85.83 


14.91 


758 


7.58 


1.05 


107.00 


3.85 


5.50 


73.88 


3.36 


17.55 


Y 


3 


2;3;4;5;6 


HD021197 


K4V 


1.50 


65 


7.36 


3.46 


65 


0.66 


0.99 


-3.50 


4.96 


-2.19 


16.25 


5.33 


2.40 


N 


3 


2;11;15 


HD022001 


F5V 


2.80 


234 


26.44 


18.02 


234 


2.34 


1.00 


28.27 


2.47 


0.74 


32.46 


1.74 


8.30 


Y 


2 


2;22 


HD022484 


F8V 


6.70 


505 


60.66 


25.91 


535 


5.36 


1.06 


108.10 


5.19 


9.14 


76.48 


2.58 


19.60 


Y 


3 


4;6;9;12 


HD022496 


K5V 


1.30 


57 


6.57 


3.11 


58 


0.59 


1.01 


6.42 


4.84 


-0.03 


0.38 


4.02 


-0.68 


N 


2 


2;3 


HD023281 


A5m 


0.39 


66 


7.54 


3.61 


97 


0.98 


1.46 


34.44 


2.63 


10.23 


14.88 


6.37 


1.77 


Y 


2 


7 


HD023754 


F5IV-V 


4.00 


389 


44.47 


23.16 


393 


3.93 


1.01 


48.81 


2.97 


1.46 


40.69 


5.62 


3.12 


N 


2 


2;3;4 
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MIPS 








PACS 




Far 






Name 


Sp. 


Age 


P 2 4 


P70 


P100 




°2*4* 


R-2A 




070 


X70 


^100 


1 

°100 


X100 


IR 


A O 

Age 


Age ref. 




type 


(Gyr) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


Exc? 


Flag 




HD027290 


F1V 


2.10 


295 


35.77 


15.60 


316 


3.16 


1.07 


220.40 


5.05 


36.56 


183.84^ 


3.64 


46.22 


Y 


1 


2 


HD029875 


F2V 


1.50 


252 


29.72 


14.43 


262 


2.62 


1.04 


-8.50 


4.46 


-8.57 


18.60 


5.88 


0.71 


N 


2 


2;16;23 


HD030652 1 " 


F6V 


1.20 


1051 


123.83 


52.97 


1093 


10.93 


1.04 


129.30 


5.53 


0.99 


91.32 


7.64 


5.02 


? 


3 


2;4;9;10;16 


HD032450 


MOV 




125 


15.03 


6.88 


133 


1.33 


- 


14.06 


3.04 


-0.32 


9.16 


4.07 


0.56 


N 






HD033111 t 


A3III 


0.39 


806 


93.08 


40.57 


822 


8.22 


1.02 


91.58 


5.56 


-0.27 


39.28 


5.16 


-0.25 


N 


2 


7 


HD033262 


F9VFe 


0.50 


311 


37.65 


16.49 


332 


3.33 


1.07 


66.32 


4.86 


5.90 


36.25 


3.03 


6.52 


Y 


3 


2;3;4;11;16 


HD033793 


sdMl.O 


1.70 


93 


11.56 


-2.22 


102 


1.03 


- 


4.28 


3.87 


-1.88 


-2.22 


5.10 




N 




2 


HD036435 


G5V 


0.50 


60 


6.87 


3.22 


61 


0.62 


1.02 


8.81 


4.40 


0.44 


6.66 


1.37 


2.51 


N 


3 


2;3;5;6;14 


HD036705 


KOV 


0.01 


99 


12.12 


5.48 


107 


1.08 


1.08 


-3.14 


6.08 


-2.51 


8.61 


0.82 


3.82 


N 


3 


2;3;13;16 


HD038393 1 " 


F6V 


2.98 


216 


44.48 


48.45 


210 


2.11 


0.97 


87.40 


2.94 


14.60 


54.72 


4.86 


1.29 


7 


B 


3;5;6 


HD038678 1 " 


A2IV-V 


0.23 


342 


38.66 


20.91 


878 


8.78 


2.57 


275.70 


3.78 


62.71 


127.63 


1.82 


58.64 


Y 


2 


7 


HD039091 


GOV 


5.80 


144 


16.30 


7.70 


144 


1.44 


1.00 


23.18 


2.98 


2.31 


8.65^ 


1.87 


0.51 


N 


3 


2;3;5;6 


HD042581 


M1/M2V 




188 


21.51 


10.05 






- 


- 






8.30 


2.47 


-0.71 


N 






HD050241 


A8VnkA6 


0.70 


641 


74.73 


32.88 


660 


6.60 


1.03 


86.12 


5.87 


1.94 


44.01 


5.51 


2.02 


N 


2 


7 


HD055892 


F3VFe 


3.90 


247 


30.53 


13.88 


270 


2.70 


1.09 


33.95 


3.11 


1.10 








N 


2 


2;23 


HD056537 


A3V 


0.55 


330 


36.94 


18.28 


326 


3.27 


0.99 


35.22 


6.60 


-0.26 


10.94 


4.83 


-1.52 


N 


2 


7 


HD056986 f 


FOIV 


4.80 


648 


74.13 


35.71 






1.05 


- 






46.52 


8.86 


1.22 


N 


2 


2;23 


HD058946 


FOV 


1.70 


350 


42.00 


18.93 


371 


3.71 


1.06 


46.39 


6.36 


0.69 


28.45 


4.76 


2.00 


N 


2 


2;16 


HD060179 


A1V 


0.25 


1826 


208.92 


112.40 






- 


- 






92.56 


4.77 


-4.16 


N 


2 


7 


HD061421* 


F5IV-V 


2.70 


13347 


1527.43 


845.64 






1.00 


- 






862.27 


3.96 


4.20 


? 


3 


2;12;15 


HD061606 


K2V 


0.55 


81 


9.06 


4.28 


80 


0.81 


0.99 


-0.77 


5.65 


-1.74 


2.63 


7.88 


-0.21 


N 


3 


2;9;10;11;15 


HD068146 


F6.5V 


3.27 


134 


15.17 


7.18 


134 


1.35 


1.00 


16.56 


2.08 


0.67 


11.63 


10.59 


0.42 


N 


1 


3 


HD071155 


AOV 


0.17 


192 


21.78 


10.60 


307 


3.08 


1.60 


211.80 


4.13 


46.01 


70.11^ 


6.59 


9.03 


Y 


2 


7 


HD071243 


F5VFe 


1.50 


418 


47.82 


22.91 






1.06 


- 






33.45 


6.06 


1.74 


N 


3 


3;16;21 


HD075632 


K5 


0.43 


117 


13.43 


6.58 






0.89 


- 






4.02 


5.69 


-0.45 


N 




2 


HD076644 


A7V 


0.75 


632 


66.52 


32.81 


651 


6.51 


1.03 


- 






37.15 


5.56 


0.78 


N 


2 


7 


HD076932 


G2VFe 


2.94 


131 


15.16 


7.10 


134 


1.35 


1.02 


14.64 


2.15 


-0.24 


11.31 


5.93 


0.71 


N 




3 


HD076943 


F4V 


1.40 


498 


58.62 


25.89 


518 


5.18 


1.04 


19.54 


26.05 


-1.50 


30.39 


6.17 


0.73 


N 


3 


2;15;16;23 


HD078045 


hA5mA5V 


0.42 


233 


27.14 


12.70 


240 


2.40 


1.03 


28.57 


4.76 


0.30 


17.45 


6.25 


0.76 


N 


2 


7 


HD078154 


F6IV 


4.90 


271 


31.00 


15.41 






0.98 


- 






8.68 


7.16 


-0.94 


N 


2 


2;12 


HD078209 


Aim 


0.80 


208 


22.86 


10.71 


202 


2.02 


0.97 


33.78 


4.09 


2.67 


19.07 


6.38 


1.31 


N 


2 


7 


HD079096 


G9V 


3.70 


98 


10.95 


5.23 


97 


0.98 


0.99 


15.72 


11.63 


0.41 


4.91 


6.32 


-0.05 


N 


2 


2;15 


HD079210 


MO.OV 


0.50 


216 


22.15 


22.98 


196 


1.96 


- 


28.22 


2.41 


2.52 


8.67 


6.39 


-2.24 


N 


1 


2 


HD079439 


A5V 


0.71 


130 


14.83 


7.18 






1.03 


- 






5.85 


6.34 


-0.21 


N 


2 


7 


HD080081 


A3V 


0.33 


255 


29.21 


14.26 






1.03 


- 






17.72 


5.32 


0.65 


N 


2 


7 


HD081997 


F6V 


1.50 


284 


32.83 


15.20 


290 


2.90 


1.02 


33.84 


3.16 


0.32 


14.52 


6.22 


-0.11 


N 


3 


1;2;4;16 


HD082328 1 " 


F7V 


5.80 


1176 


138.52 


69.12 


1223 


12.23 


1.04 


143.60 


7.26 


0.70 


66.55 


6.26 


-0.41 


N 


1 


15 


HD082885 


G8IIIv 


1.30 


243 


27.83 


13.69 






1.00 


- 






3.33 


6.56 


-1.58 


N 


3 


2;11;13;15;20 


HD084737 


G0.5Va 


9.30 


253 


28.94 


13.38 


255 


2.56 


1.01 


35.56 


3.94 


1.68 


0.24 


5.74 


-2.29 


N 


3 


9;12;15;24 


HD085376 


A5IV 


0.45 


89 


10.14 


4.99 






1.02 


- 






0.31 


6.08 


-0.77 


N 


2 


7 


HD087696 


A7V 


0.75 


175 


19.76 


9.26 


213 


2.14 


1.22 


37.46 


5.69 


3.11 


23.34 


5.77 


2.44 


Y 


2 


7 


HD088955 


A2Va 


0.41 


218 


26.17 


12.27 


231 


2.31 


1.06 


51.08 


5.89 


4.23 


22.25 


6.01 


1.66 


? 


2 


7 


HD089021 


A2IV 


0.38 


316 


36.14 


18.53 






0.95 


- 






11.63 


6.27 


-1.10 


N 


2 


7 


HD089125 


F8Vbw 


4.81 


116 


12.89 


6.39 


114 


1.15 


0.98 


22.76 


12.66 


0.78 


8.34 


6.30 


0.31 


N 


2 


4;15 


HD089269 


G5 


5.90 


70 


-193.20 


3.72 


69 


0.70 


0.99 


-193.20 


42.60 




7.55 


5.98 


0.64 


N 


3 


9;12;15;24 


HD089449 


F6IV 


3.10 


244 


27.39 


12.96 


242 


2.42 


0.99 


29.10 


3.36 


0.51 


2.24 


6.02 


-1.78 


N 


2 


1;15 


HD090132 


A8V 


0.07 


92 


10.48 


5.07 






1.05 








9.52 


6.01 


0.74 


N 


2 


7 


HD091312 


A7IV 


0.42 


144 


16.51 


8.04 






1.02 








-38.93 


11.57 


-4.06 


N 


2 


7 


HD095418 1 " 


A1V 


0.31 


837 


94.97 


50.07 


1071 


10.71 


1.28 


456.70 


9.59 


37.72 


172.94^ 


5.42 


22.67 


Y 


2 


7 


HD095608 


Aim 


0.32 


135 


15.32 


7.51 


219 


2.20 


1.62 


26.11 


7.76 


1.39 


6.92 


5.90 


-0.10 


N 


2 


7 


HD095735 


M2.0V 


2.00 


434 


56.71 


23.56 


501 


5.01 




48.23 


4.20 


-2.02 


26.11 


6.08 


0.42 


N 


1 


2;14 


HD097584 


K5 


0.80 


62 


7.12 


3.44 






0.98 








3.32 


6.21 


-0.02 


N 


2 


2;23 



to 

4^ 



o 

ffl 

H 

> 

r 
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MIPS 



PACS 



Name 


Sp. 


Age 


P 2 4 


P 70 


Pioo 


F 24 


°"2*4* 


^24 


^70 


C70 


X70 


^100 


^lOO 


Xioo 


IR 


Age 


Age ref. 




type 


(Gyr) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


Exc? 


Flag 




HD097603' 


A4V 


0.69 


902 


103.23 


40.32 


911 


9.12 


1.01 


98.49 


5.27 


-0.90 


37.21 


6.62 


-0.47 


N 


2 


7 


HD098231^ 


GOV 


0.35 


1012 


112.39 


52.85 


992 


9.92 


0.98 


109.40 


8.09 


-0.37 


48.34 


5.94 


-0.76 


N 


1 


2;5 


HD098712 


K5V 


0.25 


61 


7.02 


3.35 






1.02 








-1.86 


6.51 


-0.80 


N 


3 


3;13;14 


HD0992 1 1 


A7V(n) 


0.57 


265 


30.35 


14.74 






1.04 








12.26 


6.05 


-0.41 


N 


2 


7 


T TF\r\f\f\ A A 1 

HD099491 


K0IV 


4.10 


89 


10.40 


4.96 


92 


0.93 


1.03 


5.90 


7.38 


-0.61 


6.84 


5.87 


0.32 


N 


3 


2;9;12 


tttni aai on 

HD100180 


GOV 


3.80 


76 


8.79 


4.33 


78 


0.79 


1.02 


5.21 


6.06 


-0.59 


4.21 


6.22 


-0.02 


N 


3 


r\ r\ i r\ 1 O 10 1 A 1C 

2;9;l0;l2;l3;l4;l5 


T TTA 1 A1 1 

HD101177 


GOV 


5.23 


87 


9.94 


7.13 














-3.11 


5.99 


-1.71 


N 


3 


6;9;l2;l5 


T I r\ 1 A 1 C A 1 

HD101501 


G8V 


0.90 


267 


29.64 


9.15 


262 


2.62 


0.98 


29.01 


4.51 


-0.14 


17.15 


6.72 


1.19 


N 


3 


i./r.n.in.n.n.i a .1 c 

2;6;9;10;12;13;14;15 


T I r\ 1 A 1 C O 1 

HD101581 


K4.5Vk 


2.60 


67 


7.50 


3.62 


66 


0.67 


0.99 


6.97 


5.33 


-0.10 


-0.01 


6.15 


-0.59 


N 


2 


3;17 


HD102124 


A4V 


0.48 


114 


13.08 


6.53 






0.99 








8.72 


5.91 


0.37 


N 


2 


7 


HDl 023 65 


G2V 


6.00 


353 


41.16 


18.70 


363 


3.64 


1.03 


47.80 


4.10 


1.62 


5.92 


7.22 


-1.77 


N 


3 


3;4;5;6;12 


HDl 02647 


A3Va 


0.10 


1202 


136.69 


81.12 


1647 


16.47 


1.37 


743.00 


4.66 


130.11 








Y 


2 


7 


TTTM AOOTAT 

HDl 02870 ' 


F9V 


4.40 


881 


101.82 


46.29 


899 


8.99 


1.02 


131.20 


7.67 


3.83 


53.25 


6.69 


1.04 


N 


3 


O - A - A- 1 A- 1 O 

2;4;9;10;12 


TTTM AOOOTl 

HD 103287 1 


AOVe 


0.40 


792 


89.73 


41.73 


792 


7.92 


1.00 


95.11 


3.87 


1.39 


30.72 


5.37 


-2.05 


N 


2 


7 


T I T~\ 1 A A C 1 O 

HD104513 


A7m 


0.15 


103 


11.78 


5.82 






1.00 








-18.14 


6.62 


-3.62 


N 


2 


7 


T I TA 1 AC A CO 

HD 105452 


F1V 


1.00 


398 


45.52 


21.21 


402 


4.02 


1.01 


48.37 


4.32 


0.66 


21.65 


6.35 


0.07 


N 


3 


2;12;16 


T I r\ 1 A/C C 1 £ 

HD106516 


F9VFe 


1.00 


80 


9.12 


4.39 


80 


0.81 


1.00 


6.40 


6.18 


-0.44 


10.66 


6.21 


1.01 


N 


3 


1;3;15 


HD10o5yl 


A3V 


0.49 


426 


48.71 


22.55 


430 


4.30 


1.01 


54.18 


4.21 


1.30 


21.37 


5.92 


-0.20 


N 


2 


7 


HD 108767 


A0IVkB9 


0.26 


427 


48.37 


25.03 


427 


4.27 


1.00 


45.36 


4.86 


-0.62 


11.27 


6.31 


-2.18 


N 


2 


7 


TTTM AO AC A 

HD 108954 


F9V 


4.10 


83 


0.00 


4.42 


83 


0.83 


1.00 


0.00 


4.44 




11.31 


6.04 


1.14 


N 


2 


2;15 


TTTM (\(\(\0 C 

HD 109085 


F2V 


2.40 


377 


42.82 


19.53 


599 


6.00 


1.59 


259.10 


4.09 


52.88 


_§ 






Y 


1 


2 


HD109358^ 


GOV 


4.90 


541 


63.12 


35.64 


557 


5.57 


1.03 


60.03 


5.24 


-0.59 


31.22 


6.60 


-0.67 


N 


3 


2;9;10;12;15;20 


TTTM AACO^T 

HD 1095 36 


A7V 


0.81 


105 


12.23 


5.64 


108 


1.08 


1.03 


8.00 


3.04 


-1.39 


-4.92 


6.07 


-1.74 


N 


2 


7 


TTTM AHTO^ 

HD 1097 87 


A2V 


0.31 


235 


27.97 


12.77 


247 


2.47 


1.05 


22.73 


4.56 


-1.15 


2.31 


6.30 


-1.66 


N 


2 


7 


HDl 10304^ 


A1IV+ 


0.45 


979 


116.46 


65.46 


1028 


10.28 


1.05 


110.90 


4.28 


-1.30 


51.36 


6.13 


-2.30 


N 


2 


7 


T 1 1 1 AO 1 c 

HD1 10315 


K4.5V 


6.60 


69 


7.67 


3.72 


68 


0.69 


0.98 


19.92 


6.55 


1.87 


-3.58 


6.46 


-1.13 


N 


3 


12; 15 


TTTM 1 AO TAT 

HD1 10379 ' 


F0V 


1.18 


1359 


134.70 


67.07 


1373 


13.73 


1.01 








60.56 


6.03 


-1.08 


N 


2 


2;16 


T TF\ 1 1 A /I 1 1 

HD1 10411 


A0V 


0.50 


95 


10.70 


5.11 


147 


1.48 


1.55 


239.90 


4.96 


46.21 


140.5 


5.54 


24.44 


Y 


2 


7 


t i r\ 1 1 1 o 1 

HD111631 


K7 


0.60 


89 


9.95 


4.53 


88 


0.89 


0.99 


8.77 


6.23 


-0.19 


16.79 


5.84 


2.10 


N 


1 


2 


HD1 12758 


G9V 


8.50 


45 


5.15 


2.41 


45 


0.47 


1.01 


-5.42 


5.18 


-2.04 


2.09 


8.02 


-0.04 


N 


2 


3;15 


T I P\ 1 1 AO mo 

HD1 14378 


F5V 


0.50 


387 


43.42 


22.16 


383 


3.84 


0.99 


51.60 


5.11 


1.60 


22.48 


6.36 


0.05 


N 


3 


2;4;14;15;16 


HD1 14710 


GOV 


4.00 


512 


58.55 


25.93 


517 


5.17 


1.01 


50.47 


5.39 


-1.50 


23.79 


6.29 


-0.34 


N 


3 


1 . A. 1 A. 1 1 .10.10 

1;2;9;10;11;12;13 


HDl 15617 


G7V 


5.02 


449 


51.88 


24.86 


458 


4.58 


1.02 


156.00 


8.27 


12.59 


127.57^ 


6.15 


16.70 


Y 


3 


0/ICA1O10 1/1 

3;4;5;9;12;13;14 


HDl 15892^ 


A3mA3va 


0.26 


623 


79.08 


33.52 


698 


6.98 


1.12 


97.14 


4.55 


3.97 


29.26 


5.75 


-0.74 




2 


7 


T TTA 1 1 /I A O 

HDl 16442 


G5 


6.74 


65 


7.39 


3.54 






1.04 








5.34 


6.00 


0.30 


N 


3 


9;12;15 


HDl 16656^ 


A2V 


0.37 


872 


100.77 


82.10 


890 


8.90 


1.02 


110.20 


4.67 


2.02 


43.86 


6.29 


-6.08 


N 


2 


7 


HDl 17043 


G6V 


6.90 


82 


9.44 


4.52 






1.04 








3.59 


6.67 


-0.14 


N 


3 


12 


HDl 18098 


A3V 


0.49 


423 


48.41 


23.80 


427 


4.27 


1.01 


43.30 


5.81 


-0.88 


19.57 


5.71 


-0.74 


N 


2 


7 


HDl 18926 


K5 




31 


3.54 


1.71 






0.98 








-1.75 


6.29 


-0.55 


N 






t i r\ 1 1 a^ c 

HDl 19756 


F2V 


1.10 


326 


39.48 


21.20 


345 


3.46 


1.06 








13.50 


6.58 


-1.17 


N 


3 


2;4;16 


TTTM 1AOCA 

HDl 19850 


A It A A A 7 

M4.0V 




150 


18.65 


8.19 


165 


1.65 




18.42 


3.26 


-0.07 


20.19 


6.06 


1.98 


N 






HDl 20036 


K6+K7V 


0.80 


45 


5.26 


2.41 


46 


0.48 


1.02 


10.27 


6.11 


0.82 


-3.51 


6.50 


-0.91 


N 


1 


2;3 


HDl 20467 


K5.5Vk 


4.35 


65 


7.29 


3.40 


64 


0.65 


0.98 


-3.69 


7.04 


-1.56 


-8.04 


5.75 


-1.99 


N 


3 


3;12 


nL'lZ'tjoU 






oU 


O. ly 




oU 


U.ol 


l.Ul 








LJ.J 1 


0. 1U 


1 AQ 


XT 
In 


j 


J,D,V 


HD125161 


A7V 


0.04 


139 


16.35 


7.62 


144 


1.45 


1.04 


18.50 


3.52 


0.61 


15.42 


6.19 


1.26 


N 


2 


1 


HD125162 


AOp 


0.29 


192 


21.79 


5.73 


282 


2.83 


1.47 


378.30 


6.66 


53.53 


240.49^ 


5.55 


42.30 


Y 


2 


1 


HD126660 1 " 


F7V 


0.50 


546 


64.32 


35.29 


568 


5.68 


1.04 


70.77 


5.16 


1.25 


26.01 


5.91 


-1.57 


N 


3 


2;4;16;23 


HD128167 


F2V 


1.70 


280 


32.00 


15.62 






1.02 








29.62 


3.89 


3.60 


? 


2 


2;12 


HDl 29502 


F2V 


1.80 


520 


61.20 


27.28 


540 


5.41 


1.04 


55.87 


3.65 


-1.46 


25.49 


5.98 


-0.30 


N 


1 


2;23 


HD130109 


A0V 


0.29 


246 


27.32 


13.22 


241 


2.41 


0.98 


-12.06 


4.75 


-8.29 


-0.26 


5.86 


-2.30 


N 


2 


7 


HD130841 t 


A5mA4IV 


0.65 


154 


86.14 


45.19 


143 


1.44 


0.93 








34.97 


6.55 


-1.56 


N 


2 


7 


HD131156 


G8V 


0.28 


489 


55.95 


27.26 






1.01 








40.68 


5.99 


2.24 


N 


3 


2;5;9;10;12;13;14;15 
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Table 2 — Continued 



to 



MIPS 



PACS 



Name 


Sp. 


Age 


P 2 4 


P 70 


Pioo 


F 24 


°"2*4* 


^24 


^70 


C70 


X70 


^100 


^lOO 


Xioo 


IR 


Age 


Age ref. 




type 


(Gyr) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


Exc? 


Flag 




HD 133640 


GOVnv 


0.95 


410 


46.94 


22.06 






1.09 








34.70 


5.99 


2.11 


N 


1 


23 


HD 135204 


K0V 


8.57 


88 


10.08 


4.80 


89 


0.90 


1.01 


6.09 


5.96 


-0.67 


8.71 


5.67 


0.69 


N 


1 


15 


HD 136923 


G9V 


3.00 


54 


6.07 


2.95 


54 


0.55 


1.00 


11.70 


4.54 


1.24 








N 


3 


9;12;15;20 


HD137107 


GOV 


2.60 


248 


28.42 


13.50 


251 


2.51 


1.01 


42.00 


5.78 


2.35 


11.08 


5.91 


-0.41 


N 


1 


2 


HD 137763 


G9V 


6.67 


80 


8.81 


4.31 


80 


0.81 


1.00 








3.72 


5.91 


-0.10 


N 


3 


9;10;15 


HD 137898 


A8IV 


0.07 


104 


11.65 


5.53 


103 


1.04 


0.99 


15.64 


3.87 


1.03 


-7.07 


6.63 


-1.90 


N 


2 


7 


HD 137909 


FOp 


0.81 


366 


38.21 


19.25 


358 


3.59 


0.98 








16.06 


6.25 


-0.51 


N 


2 


7 


HD 139006 1 


A0V 


0.27 


983 


111.10 


59.83 


1298 


12.98 


1.32 


509.50 


9.82 


40.57 


211.48^ 


6.21 


24.42 


Y 


2 


7 


HD 139763 


K6Vk 


1.60 


40 


4.65 


2.16 


42 


0.44 


1.04 








6.36 


6.17 


0.68 


N 


3 


2;3;11 


HD 140436 


B9IV+ 


0.40 


215 


25.08 


11.66 


222 


2.22 


1.03 


23.11 


4.02 


-0.49 


2.69 


6.55 


-1.37 


N 


2 


7 


T i r\ 1/11 rvA/i 

HD 14 1004 


G0IV-V 


5.30 


466 


50.70 


25.71 


447 


4.48 


0.96 


52.38 


5.10 


0.33 


25.76 


5.44 


0.01 


N 


3 


2;4;5;9;10;12;13;14;15 


T TTA 1/11 

HD141272 


G8V 


0.70 


44 


5.04 


2.42 


44 


0.47 


1.00 


-17.41 


5.41 


-4.15 


3.21 


5.67 


0.14 


N 


3 


2;9;10;12;13;14;15 


HD141795 


A2m 


0.52 


280 


32.08 


15.68 






1.03 








18.79 


5.86 


0.53 


N 


2 


7 


HD 142267 


GOV 


4.80 


109 


12.48 


5.96 


110 


1.11 


1.01 


10.78 


1.72 


-0.99 


6.37 


5.92 


0.07 


N 


3 


2;9;10;12 


HD142373 


F8Ve 


6.21 


427 


49.40 


22.77 


436 


4.36 


1.02 


40.25 


5.58 


-1.64 


9.45 


6.00 


-2.22 


N 


3 


9;12;15 


HD 1428 60^ 


F6IV 


4.60 


640 


75.38 


32.05 


665 


6.66 


1.04 


73.56 


6.51 


-0.28 


27.36 


6.80 


-0.69 


N 


3 


2;4;9;l0;l2 


HD146361 


F6V+G0V 


0.01 


179 


21.34 


9.22 


188 


1.89 


1.05 


28.28 


5.22 


1.33 


7.36 


5.65 


-0.33 


N 


3 


2;15;16 


HD 147379 


M1V 




79 


9.41 


4.36 


83 


0.83 




17.10 


5.06 


1.52 


1.47 


1.13 


-2.56 


N 






HD 1475 84 


F9V 


1.70 


296 


30.73 


15.32 


288 


2.88 


0.97 








14.64 


6.20 


-0.11 


N 


3 


3;5;11;14 


HD151288 


K5 


2.60 


103 


11.58 


5.35 


102 


1.03 


0.99 


16.45 


1.65 


2.95 


0.18 


6.08 


-0.85 


N 


3 


2;l2;l5;23 


HD 154494 


A4IV 


0.40 


97 


11.14 


5.46 






1.01 








2.76 


5.75 


-0.47 


N 


2 


7 


HD 154577 


K2.5Vk 


4.83 


66 


7.27 


3.54 


66 


0.67 


1.00 








2.32 


6.41 


-0.19 


N 


2 


3;l7 


TTTA1 CCO'lH 

HD155876 


K5 




118 


13.53 


6.36 






1.05 








1.74 


2.42 


-1.91 


N 






HD156164 


A3IV 


0.35 


513 


58.65 


26.71 


518 


5.18 


1.01 


57.45 


4.81 


-0.25 


28.84 


5.76 


0.37 


N 


2 


7 


HD 159560 


A4m 


0.70 


145 


16.62 


8.15 






1.01 








13.66 


6.26 


0.88 


N 


2 


7 


tttai zrrvr\o / "> 

HD 160032 


F4V 


2.22 


236 


27.26 


12.70 


241 


2.41 


1.02 


47.44 


4.35 


4.64 


37.03 


6.72 


3.62 


Y 


1 


2 


HD 160922 


F4V 


2.50 


129 


27.64 


7.24 


244 


2.44 




32.13 


5.22 


0.86 


8.41 


5.83 


0.20 


N 


1 


2;23 


HD 162003 


F5IV-V 


4.20 


300 


34.36 


16.06 


303 


3.03 


1.01 


47.08 


4.97 


2.56 


17.60 


6.41 


0.24 


N 


2 


l;2 


T TTT\ 1 H Cf\ Af\ 

HD 165040 


A7sp 


0.80 


221 


24.48 


11.57 


216 


2.16 


0.98 


-3.83 


4.89 


-5.79 


17.61 


6.04 


1.00 


N 


2 


7 


t i r\ 1 nc 1 oc\ 

HD165189 


A5V 


0.01 


118 


13.53 


6.61 






1.01 








9.32 


6.60 


0.41 


N 


3 


16;25 


TTTA1 HC~IT~I 

HD 165777 


A4IVs 


0.55 


275 


31.49 


15.26 






1.04 








17.17 


6.38 


0.30 


N 


2 


7 


HD 165908 


F7V 


7.20 


273 


30.64 


14.61 


270 


2.70 


0.99 


98.65 


5.18 


13.13 


78.42^ 


6.64 


9.61 


Y 


3 


2;9;l0 


TTTA1 HdlO AO 

HD 166348 


MOV 


1.65 


61 


7.14 


3.34 


63 


0.64 




19.24 


4.84 


2.50 


3.61 


6.77 


0.04 


N 


1 


3 


TTTA1 /T/IOC 

HD 167425 


F9.5V 


0.90 


77 


8.83 


4.25 






1.04 








4.89 


5.79 


0.11 


N 


3 


2;3;5 


T TT\ 1 HO 1 C 1 

HD168151 


F5V 


2.50 


212 


24.20 


11.16 


214 


2.14 


1.01 


24.49 


3.21 


0.09 


6.73 


6.15 


-0.72 


N 


2 


l;2;23 


HD170153 ' 


F7V 


5.50 


997 


115.21 


60.05 


1017 


10.17 


1.02 


129.80 


5.05 


2.89 


46.56 


5.62 


-2.40 


N 


2 


2;23 


HD172555 


A7V 


0.01 


136 


15.41 


13.93 


866 


8.66 


6.37 


226.40 


5.95 


35.46 


81.52 


6.04 


11.19 


Y 


2 


7 


t i r"\ 1 nonon 

HD173739 


M3.0V 




101 


21.30 


5.70 


188 


1.88 




39.20 


3.53 


5.07 


19.50 


3.04 


4.54 


Y 






HD 176051 


GOV+klV 


3.50 


254 


27.93 


13.76 


247 


2.47 


0.97 


27.12 


5.09 


-0.16 


15.03 


6.34 


0.20 


N 


2 


2;l4 


TTTA1 'I HHO'IX 

HD 176687 1 


A2.5Va 


0.48 


796 


87.49 


46.46 


773 


7.73 


0.97 


65.60 


7.32 


-2.99 


31.00 


5.97 


-2.59 


N 


2 


7 


HD177196 


A7V 


0.59 


104 


11.93 


5.85 






1.02 








12.97 


6.14 


1.16 


N 


2 


7 


TTTA1 1 (\ 

HD 179930 


A KCW 71 

MOVk 


1.49 


49 


5.56 


2.70 














6.95 


6.25 


0.68 


N 


1 


3 


T TT\ 1 O (\ 1 1 

HD180161 


G8V 


0.60 


60 


6.77 


3.26 


60 


0.61 


0.99 


8.62 


4.41 


0.42 


15.13 


6.25 


1.90 


N 


3 


2;9;12;13;14 


M1J 1 oU / / / 


A7V 


n 1 n 

U. 1U 




14 71 


/ .uo 














^ 1 7 

0. 1 / 


A ^7 
D.J 1 


-i.yz 


XT 
In 


2 


/ 


HD181321 


G2V 


0.15 


79 


9.16 


4.26 


81 


0.83 


1.02 


0.51 


2.50 


-3.46 


1.90 


5.63 


-0.42 


N 


3 


2;3;5;6;16 


HD 184006 


A5V 


0.45 


358 


38.91 


18.77 


343 


3.44 


0.96 


38.67 


6.01 


-0.04 


13.40 


5.54 


-0.97 


N 


2 


7 


HD186219 


A4III 


0.58 


84 


9.56 


4.62 






1.03 








11.88 


6.15 


1.18 


N 


2 


7 


HD186408 


G1.5Vb 


7.54 


110 


11.11 


6.50 


98 


0.99 


0.89 


-1.13 


2.90 


-4.22 


6.86 


6.06 


0.06 


N 


3 


9;10;12;15 


HD1 87642 f 


A7V 


0.70 


4558 


521.65 


298.00 






1.04 








292.02 


3.96 


-1.51 


N 


2 


7 


HD1 88228 


AOVa 


0.25 


172 


20.25 


9.08 


179 


1.79 


1.04 


73.47 


5.90 


9.02 


42.00^ 


6.13 


5.37 


Y 


2 


7 


HD 189245 


F7V 


0.06 


120 


14.10 


6.40 


125 


1.26 


1.04 


10.59 


2.34 


-1.50 


2.48 


6.12 


-0.64 


N 


3 


2;3;4;16 


HD 190007 


K4Vk: 


1.80 


90 


10.23 


4.79 


90 


0.92 


1.00 


10.68 


3.01 


0.15 


1.28 


6.39 


-0.55 


N 


3 


2;3;12;13;14;15;20 



o 

GO 
ffl 

H 

> 

r 



Table 2 — Continued 



MIPS PACS Far 



Name 


Sp. 


Age 


P 2 4 


P 70 


Pioo 


F 24 


°"2*4* 


^24 


^70 


C70 


X70 


^100 


^lOO 


Xioo 


IR 


Age 


Age ref. 




type 


(Gyr) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


Exc? 


Flag 




HD 190422 


F9V 


0.40 


74 


8.73 


4.04 


77 


0.78 


1.04 


-0.38 


6.37 


-1.43 


5.05 


6.29 


0.16 


N 


3 


2;3;5;16 


HD191849 


MOV 


1.00 


169 


20.60 


9.32 


182 


1.83 




33.79 


2.38 


5.54 


32.05 


5.71 


3.98 


Y 


1 


2;17 


HD 1923 10 


K2+V 


6.10 


260 


29.19 


13.58 


258 


2.58 


0.99 


23.74 


3.68 


-1.48 


11.39 


6.64 


-0.33 


N 


3 


2;3;12;17 


T i r\ 1 C\ A f Af\ 

HD 194640 


G8V 


4.92 


76 


8.50 


4.12 


75 


0.76 


0.99 


7.18 


6.29 


-0.21 


3.93 


6.29 


-0.03 


N 


3 


3;5;6 


HD 196877 


K7V 


4.55 


50 


6.11 


2.82 


54 


0.56 


1.08 


3.55 


1.91 


-1.34 


-2.61 


5.72 


-0.95 


N 


1 


3 


HD 197076 


G5V 


4.75 


76 


8.47 


4.10 


75 


0.75 


0.98 


-3.96 


5.20 


-2.39 








N 


3 


6;9;12;15;20 


HD197157 


A9IV 


0.52 


191 


21.91 


10.69 






1.03 








14.35 


6.20 


0.59 


N 


2 


7 


HD 197692 


F5V 


1.00 


424 


48.01 


24.23 


424 


4.25 


1.00 


-6.95 


7.56 


-7.27 


19.17 


6.40 


-0.79 


N 


3 


2;3;4;16;23 


HD200525 


F9.5V 


1.00 


157 


17.94 


8.84 






1.00 








10.56 


5.72 


0.30 


N 


3 


2;3;5 


HD200779 


K6V 


5.55 


54 


6.20 


3.06 






0.97 








8.55 


5.78 


0.95 


N 


2 


1 1 ; 15 


T I r\0 C\f\C\ £ o 

HD200968 


K1IV 


1.25 


68 


7.51 


-2.67 


66 


0.68 


0.97 


-3.04 


5.96 


-1.77 


-2.67 


6.30 




N 


3 


3;11;12 


HD202275 


F5V+ 


4.90 


363 


41.47 


18.69 


366 


3.66 


1.01 


45.14 


4.64 


0.79 


20.13 


6.25 


0.23 


N 


1 


15 


HD202560 


M0V 


4.78 


511 


58.50 


29.63 














36.70 


8.13 


0.87 


N 


2 


17 


HD202730 


A5V(n) 


0.60 


162 


18.14 


7.68 


160 


1.61 


0.99 


24.99 


4.60 


1.49 


20.32 


6.32 


2.00 


N 


2 


7 


HD203244 


G5V 


0.39 


60 


6.90 


3.33 






1.03 








7.41 


6.09 


0.67 


N 


3 


2;3;5;6 


HD203608 


F9VFe 


0.57 


493 


57.48 


25.54 


507 


5.07 


1.03 


51.86 


6.11 


-0.92 


26.27 


6.05 


0.12 


N 


1 


3 


HD204961 


Ml. 5 




155 


19.16 


8.46 


169 


1.69 




22.62 


4.95 


0.70 


27.39 


5.86 


3.23 


N 






HD206826 


F6V 


2.70 


335 


38.30 


18.07 


338 


3.38 


1.01 


39.36 


4.80 


0.22 


26.97 


5.97 


1.49 


N 


2 


2;15 


HD207098 ' 




0.01 


1012 


110.50 


54.60 


1002 


10.02 


0.99 








57.93 


6.06 


0.55 


N 


3 


26 


HD2 10027^ 


F5V 


5.20 


643 


75.00 


43.68 


662 


6.62 


1.03 


71.16 


5.26 


-0.73 


28.63 


6.02 


-2.50 


N 


2 


2;23 


HD2 10049 


A1.5IVn 


0.39 


124 


14.23 


7.01 






1.03 








12.74 


6.16 


0.93 


N 


2 


7 


HD210418 


AlVa 


0.50 


341 


39.43 


18.32 


348 


3.48 


1.02 


46.92 


4.83 


1.55 


10.89 


6.88 


-1.08 


N 


2 


7 


I I r\0 1 1 C\'~IC\ 

HD2 11970 


K7Vk 


1.70 


44 


5.01 


2.40 


44 


0.45 


1.00 


11.30 


5.77 


1.09 


1.91 


6.11 


-0.08 


N 


1 


3 


HD212330 


G2IV-V 


7.90 


242 


27.90 


12.92 


247 


2.47 


1.02 


-5.08 


4.34 


-7.60 


10.15 


6.29 


-0.44 


N 


3 


3;5;6;17;24 


HD2 12698 


G2V 


0.35 


164 


18.79 


9.13 






1.01 








12.47 


6.30 


0.53 


N 


3 


2;3;11 


HD212728 


A4V 


0.06 


66 


7.55 


3.69 






1.02 








3.26 


6.12 


-0.07 


N 


2 


7 


HD213398 


A0V 


0.18 


139 


15.71 


7.74 


171 


1.71 


1.23 


63.48 


3.52 


13.57 


31.48 


6.01 


3.95 


Y 


2 


7 


HD2 13845 


F7V 


0.90 


157 


18.32 


8.37 


162 


1.63 


1.03 


-1.42 


2.96 


-6.67 


15.05 


6.02 


1.11 


N 


3 


2;3;4;16 


I i r\0 1 A 'I Af\ 

HD2 14749 


K4.5Vk 


0.60 


72 


8.04 


3.74 


71 


0.72 


0.98 


8.22 


5.88 


0.03 


2.74 


6.23 


-0.16 


N 


3 


2;3;12 


HD2 14953 


F9.5V 


5.17 


105 


12.02 


5.77 






1.06 








2.49 


5.65 


-0.58 


N 


3 


3;5;6 


T I r\o 1 CH AO 

HD2 15648 


F7V 


5.00 


485 


54.95 


25.95 


485 


4.85 


1.00 


54.65 


3.69 


-0.08 


13.34 


6.93 


-1.82 


N 


3 


9;12;24 


HD215789 


A2IVnSB 


0.60 


407 


42.36 


24.80 


374 


3.74 


0.92 


41.94 


4.24 


-0.10 


28.28 


6.00 


0.58 


N 


2 


7 


HD216133 


M0.5V 




31 


3.51 


1.68 














3.86 


5.74 


0.38 


N 


" 




HD2 16803 


K4V 


0.25 


230 


25.00 


12.46 


221 


2.21 


0.96 


26.37 


3.11 


0.44 


4.59 


6.25 


-1.26 


N 


3 


2;3;5;6;12;13;14;1 


HD2 16899 


M2.0V 




133 


16.05 


7.20 


142 


1.43 




18.27 


3.37 


0.66 


16.03 


5.81 


1.52 


N 






HD217107 


G8IV 


8.48 


111 


12.51 


6.02 


110 


1.12 


0.99 


6.32 


5.63 


-1.10 


1.20 


6.34 


-0.76 


N 


3 


9;10;12;15 


HD2 17987 


M2V 




404 


51.46 


28.28 


454 


4.55 




54.30 


4.65 


0.61 


28.28 


6.06 




N 






t i r\o i oci 1 

HD218511 


K6V 


0.90 


56 


6.26 


2.92 


55 


0.56 


0.99 


0.30 


5.09 


-1.17 


5.35^ 


6.24 


0.39 


N 


2 


2;3 


HUziyj /I 


F4V 


4.70 


503 


57.51 


28.33 


508 


5.08 


1.01 


54.46 


4.92 


-0.62 


25.92 


5.74 


-0.42 


N 


1 


2 


HD222335 


G9.5V 


3.36 


57 


6.43 


3.10 


57 


0.58 


0.99 


0.86 


5.51 


-1.01 


-0.43 


6.19 


-0.57 


N 


3 


3;5;6;12;17 


HD222345 


A7IV 


0.60 


127 


14.58 


7.15 






1.02 








-2.22 


6.55 


-1.43 


N 


2 


7 


HD222368 


F7V 


5.20 


536 


60.08 


26.95 


531 


5.31 


0.99 


70.05 


5.39 


1.85 








N 


2 


2;12 


HD222603 


A7V 


0.70 


167 


19.12 


4.76 






1.01 








20.46 


5.75 


2.73 


N 


2 


7 


HD223352 


A0V 


0.22 


108 


12.23 


6.10 


160 


1.61 


1.48 


54.80 


6.62 


6.43 


7.33 


5.87 


0.21 


? 


2 


7 


HD224953 


MOV 




40 


4.53 


2.19 














14.19 


6.22 


1.93 


N 






HD234078 


K5 


0.85 


50 


5.53 


2.66 


49 


0.49 


0.97 


10.88 


5.10 


1.05 


2.00 


5.49 


-0.12 


N 




15 


HD265866 


M4.0V 




78 


8.96 


4.17 


88 


0.88 










3.56 


7.63 


-0.08 


N 






DUNES survey 


HIP000171 


G5Vb 


4.00 


209 


24.40 


11.18 


215 


2.15 


1.03 


27.50 


4.92 


0.63 


10.49 


3.47 


-0.20 


N 


3 


1;2;15;23 


HIP000544 


K0V 


0.24 


137 


15.17 


7.19 


158 


1.53 


1.15 


106.00 


3.46 


26.25 


53.31 


2.71 


17.02 


Y 


3 


1;2;6;15;23 


HIP000910 


F8VFe 


3.00 


264 


30.30 


13.91 


267 


3.89 


1.01 


37.40 


3.66 


1.94 


17.90 


5.54 


0.72 


N 


3 


3;4;23;24 


HIP002941 


G8V 


5.10 


207 


23.17 


10.99 


205 


2.05 


0.99 


19.20 


5.03 


-0.79 


7.86 


4.41 


-0.71 


N 


3 


1;3;5;9;24;27 
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Sp. 


Age 


P 2 4 


P 70 


Pioo 


F 24 


°24* 




^70 


C70 


X70 


^100 


^100 


Xioo 


IR 


Age 


Age ref. 




type 


(Gyr) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


Exc? 


Flag 




HIP003093 


KOV 


6.50 


196 


22.26 


10.42 


196 


1.97 


1.00 


11.80 


4.61 


-2.27 


7.68 


3.97 


-0.69 


N 


3 


1 O C 1 A 10 1 C 

I;2;5;l0;l3;l5 


H1P003497 


G6VFe 


5.70 


75 


8.60 


4.02 


76 


0.76 


1.01 


6.69 


4.67 


-0.41 


6.84^ 


2.85 


0.99 


N 


3 


2;3;5;12;15;23;24 


HIP003821 


G3V 


5.40 


1198 


130.21 


105.97 


1150 


11.50 


0.96 


125.00 


5.85 


-0.89 


51.44 


3.84 


-14.20 


N 


3 


1;2;5;10;12;24 


HIP003909 


F7IV-V 


4.00 


197 


22.22 


10.35 


197 


1.98 


1.00 


25.50 


2.76 


1.19 


16.70 


4.38 


1.45 


N 


3 


2;4;12;15;23;24 


TTTDhA/l 1 /to 

H1P004148 


K2.5Vk 


1.55 


84 


9.20 


4.21 


80 


0.81 


0.95 


34.70 


5.10 


5.00 


18.61 ^ 


2.77 


5.20 


Y 


3 


2;3;17 


H1P007513 


F9V 


4.00 


543 


61.52 


27.42 


543 


5.43 


1.00 


56.30 


5.07 


-1.03 


30.87 


4.86 


0.71 


N 


3 


1;2;5;10;12;23;24 


HIP007978 


F9V 


1.90 


158 


26.16 


8.04 


196 


1.96 


1.24 


1 f \ a f\ r\r\ 

1040.00 


5.80 


174.80 


§ 






Y 


3 


3;4;17 


HIP008768 


MOV 


0.60 


72 


14.10 


3.66 


71 


0.71 


0.98 


14.10 


2.78 




3.75^ 


2.89 


0.03 


N 


2 


2;13 


HIP010138 


G9V 


2.20 


169 


18.81 


8.78 


166 


1.67 


0.98 


3.44 


5.61 


-2.74 


4.92 


3.98 


-0.97 


N 


3 


2;3;4;5;17 


HIPO 10798 


G8V 


4.50 


112 


12.70 


5.99 


112 


1.13 


1.00 


10.70 


2.32 


-0.86 


3.96 


2.47 


-0.82 


N 


3 


3;5;9;10;12 


HIP0 11452 


M1V 




73 


-4.90 


3.67 


69 


0.71 


0.95 


-4.90 


7.01 




11.00^ 


2.41 


3.04 


? 


" 




HIPO 11964 


K7V 


0.00 


107 


11.62 


5.18 


103 


1.04 


0.96 


1.16 


3.19 


-3.28 


8.40 


2.75 


1.17 


N 


3 


2;3;16 


HIPO 12777 


F7V 


6.00 


492 


57.44 


27.06 


507 


5.07 


1.03 


52.80 


4.64 


-1.00 


24.03 


5.61 


-0.54 


N 


3 


r\ c OO O A 

2;5;23;24 


HIPO 13402 


K1V 


0.19 


187 


22.13 


9.92 


194 


1.95 


1.04 


64.70 


6.17 


6.90 


48.50 


2.65 


14.56 


Y 


3 


1;2;3;5;12;13;17;23 


HIPO 14954 


F8V 


4.50 


242 


27.45 


13.42 


244 


2.45 


1.01 


45.50 


3.54 


5.10 


32.26 


3.12 


6.04 


Y 


3 


4;9;10;12;15;24 


HIPO 15330 


G4V 


2.00 


193 


22.15 


10.21 


195 


1.95 


1.01 


31.30 


3.66 


2.50 


-5.09 


4.07 


-3.76 


N 


3 


2;3;5;6;11;17 


HIPO 15 371 


GOV 


4.00 


230 


26.90 


12.19 


237 


2.37 


1.03 


41.70 


3.29 


4.50 


37.59^ 


3.01 


8.44 


Y 


3 


2;5;15;17;28 


HIPO 15799 


KOV 


2.30 


95 


10.84 


5.29 


94 


2.65 


0.99 








5.40^ 


2.67 


0.04 


N 


2 


3;5 


HIPO 16 134 


K7V 


0.50 


77 


8.32 


3.91 


74 


0.74 


0.96 


8.32 


4.92 




4.33 


3.01 


0.14 


N 


2 


2;3 


HIPO 17420 


K2V 


2.00 


88 


10.29 


4.45 


83 


0.84 


0.94 


24.30 


5.39 


2.60 


19.35^ 


2.79 


5.34 


Y 


1 


2;3 


HIPO 17439 


K2V 


0.60 


78 


8.75 


4.11 


77 


0.78 


0.99 


89.10 


4.32 


18.60 


75.02^ 


2.74 


25.88 


Y 


3 


2;3;17 


HIPO 19849 


K0.5V 


5.50 


761 


89.62 


39.63 


791 


7.91 


1.04 


86.60 


2.54 


-1.19 


35.70 


3.61 


-1.09 


N 


3 


1;2;3;5;9;10;12;13;15;23;27 


H1P019884 


K4.5Vk 


4.50 


77 


8.78 


4.30 


77 


0.78 


1.00 


13.20 


4.46 


0.99 


2.44 


2.95 


-0.63 


N 


2 


3;17 


H1P022263 


G1.5 


0.70 


186 


23.28 


9.95 


193 


1.94 


1.04 


120.00 


3.72 


26.00 


69.95 


2.71 


22.14 


Y 


3 


1;2;3;4;5;9;10;12;13;23;27 


HIP023311 


K3V 


5.50 


244 


26.02 


12.00 


229 


2.30 


0.94 


25.20 


2.75 


-0.30 


12.28 


5.52 


0.05 


N 


3 


1;2;10;15 


HIP025110 


F6V 


2.00 


194 


21.99 


10.15 


194 


1.94 


1.00 


-5.10 


5.13 


-5.28 


-57.21 


5.93 


-11.36 


N 


2 


1;2;4;15;16;24 


HIP027887 


K2.5V 


5.50 


94 


10.02 


4.75 


88 


0.89 


0.94 


15.00 


4.70 


1.06 


11.53^ 


2.68 


2.53 


N 


2 


3;17 


HIP028103 


F2V 


2.50 


506 


62.03 


26.58 


567 


5.67 


1.12 


96.00 


3.86 


8.80 


39.52 


6.19 


2.09 


Y 


1 


2;23;24 


HIP028442 


K6.5V 


5.30 


75 


8.62 


4.14 


75 


6.42 


1.00 








0.28^ 


2.82 


-1.37 


N 


1 


3 


HIP02927 1 


G7V 


5.00 


329 


36.45 


17.39 


322 


3.22 


0.98 


39.10 


5.41 


0.49 


4.36 


5.50 


-2.37 


N 


3 


2;3;5 


HIP029568 


G5V 


0.35 


95 


11.02 


5.12 


97 


0.98 


1.02 


14.90 


2.44 


1.59 


9.08 


2.89 


1.37 


N 


3 


2;3;9;10;12 


HIP032439 


F7V 


4.50 


158 


17.33 


8.22 


153 


1.54 


0.97 


20.50 


2.22 


1.43 


6.00 


4.26 


-0.52 


N 


3 


1;2;6;15;24 


HIP032480 


GOV 


5.00 


192 


33.00 


10.37 


200 


2.01 


1.04 


297.00 


3.30 


80.00 


192.13^ 


2.66 


68.33 


Y 


3 


1;2;4;5;12;15 


HIP033277 


GOV 


5.80 


140 


15.86 


7.24 


140 


1.41 


1.00 


8.91 


4.04 


-1.72 


8.26^ 


3.10 


0.33 


N 


3 


1;6;9;12;15;24 


HIP034017 


G4V 


6.40 


116 


13.40 


6.22 


118 


1.20 


1.02 


13.00 


3.66 


-0.11 


11.17 


3.32 


1.49 


N 


3 


1;4;6;9;12;15 


HIP034065 


GOV 


6.10 


172 


19.23 


9.53 


170 


1.71 


0.99 


23.20 


2.53 


1.57 


2.70 


4.04 


-1.69 


N 


3 


5;15;24;28 


HIP035136 


GOV 


5.70 


167 


19.29 


8.87 


170 


1.71 


1.02 


27.80 


3.96 


2.15 


4.54 


4.29 


-1.01 


N 


3 


9;12;15;24 


HIP036439 


F6V 


5.00 


155 


17.28 


8.18 


153 


1.54 


0.99 


17.00 


3.52 


-0.08 


11.03 


3.13 


0.91 


N 


2 


1;24 


HIP038382 


GOV 


5.50 


251 


28.16 


12.99 


248 


2.48 


0.99 


17.10 


5.50 


-2.01 


11.02 


5.48 


-0.36 


N 


3 


3;24 


HIP038784 


G8V 


4.10 


84 


9.41 


4.51 


83 


0.85 


0.99 


10.20 


1.97 


0.40 


5.40^ 


2.98 


0.30 


N 


3 


1;2;6;9;10;12 


HIP040693 


G8+V 


6.00 


159 


17.41 


8.99 


235 


2.36 


1.48 


15.90 


2.12 


-0.71 


7.88 


3.69 


-0.30 


N 


3 


2;3;4;9;10;12 


HIP040843 


F6V 


4.60 


203 


22.52 


10.77 


199 


2.00 


0.98 


33.80 


5.37 


2.10 


29.71 


4.53 


4.18 


Y 


3 


1;9;12;15;20 


HIP042430 


G5IV 


7.00 


281 


34.32 


15.70 


303 


3.03 


1.08 


33.80 


5.16 


-0.10 


16.96 


5.73 


0.22 


N 


3 


2;3;4;24 


rllrU4Z4o5 
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j 
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HIP043587 


G8V 


7.10 


181 


20.06 


9.22 


177 


1.78 


0.98 


19.80 


3.19 


-0.08 


10.64 


3.56 


0.40 


N 


3 


6;9;10;12;15 


HIP043726 


G3V 


2.20 


118 


13.82 


6.23 


122 


1.23 


1.03 


32.90 


3.18 


6.00 


14.84 


2.86 


3.01 


Y 


3 


1;2;3;5;12;13 


HIP044897 


F9V 


0.80 


109 


12.35 


5.77 


109 


1.10 


1.00 


17.30 


3.28 


1.51 


11.37^ 


2.87 


1.95 


N 


3 


1;2;10;12;13;15;23 


HIP045333 


GOV 


7.20 


245 


27.82 


12.77 


245 


2.45 


1.00 


24.60 


4.95 


-0.65 


13.19 


4.69 


0.09 


N 


3 


4;15;23;24 


HIP045617 


K3V 


2.00 


94 


9.96 


4.71 


88 


0.89 


0.94 


6.96 


11.10 


-0.27 


2.39^ 


2.83 


-0.82 


N 


2 


2;15 


HIP046580 


K3V 


0.45 


93 


10.15 


4.77 


90 


0.91 


0.96 


-2.19 


6.36 


-1.94 


10.45^ 


2.64 


2.15 


N 


3 


1;2;6;15 


HIP047592 


F8V 


4.50 


256 


28.98 


13.48 


256 


2.57 


1.00 


0.97 


7.45 


-3.76 


-56.61^ 


6.49 


-10.80 


N 


3 


3;4;12;17;24 


HIP049081 


G3Va 


8.10 


215 


23.30 


11.18 


211 


2.11 


0.98 


-330.00 


58.30 


-6.06 


7.05 


5.58 


-0.74 


N 


3 


1;9;12;15;24 
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Table 2 — Continued 



















MIPS 








PACS 




Far 






Name 


Sp. 


Age 


P 2 4 


P 70 


Pioo 


^24 


^2*4* 


^24 




C70 


X70 


^100 


j 

^lOO 


X100 


IR 


Age° 


Age ref. 




type 


(Gyr) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


Exc? 


Flag 




HIP049908 


K8V 


1.20 


463 


49.23 


23.03 


435 


4.35 


0.94 


43.20 


4.22 


-1.43 


17.26 


2.96 


-1.95 


N 


3 


1;2;15 


HIP051459 


F8V 


3.10 


287 


31.85 


14.68 


281 


2.81 


0.98 


34.00 


3.84 


0.56 


27.94 


6.14 


2.16 


N 


3 


1;2;4;9;10;15 


HIP051502 


F2V 


1.50 


143 


17.05 


7.65 


149 


1.49 


1.04 


-0.31 


2.17 


-8.00 


46.65 


2.82 


13.83 


Y 


2 


16;23;24 


HIP053721 


G1V 


6.50 


278 


30.58 


14.68 


270 


2.70 


0.97 


33.00 


4.24 


0.57 


5.78 


5.78 


-1.54 


N 


3 


1;9;10;12 


HIP054646 


K8V 


1.20 


81 


12.20 


4.17 


76 


0.77 


0.94 


12.20 


4.88 




8.51^ 


2.95 


1.47 


N 




2 


HIP056452 


KOV 


4.70 


189 


21.19 


9.92 


187 


1.88 


0.99 


24.90 


2.41 


1.54 


8.93 


3.95 


-0.25 


N 


3 


2;3;4;5;9;10;12 


HIP057507 


G6V 


5.10 


87 


9.81 


4.62 


87 


0.88 


1.00 


9.96 


2.16 


0.07 


6.54^ 


2.87 


0.67 


N 


2 


3;5 


HIP057939 


G8Vp 


4.50 


133 


14.69 


6.96 


130 


1.31 


0.98 


10.10 


1.89 


-2.43 








N 


3 


1;2;9;10;12;15 


HIP058345 


K4+V 


2.50 


138 


16.29 


7.55 


143 


1.43 


1.04 


10.20 


4.91 


-1.24 


7.15 


3.09 


-0.13 


N 


1 


2;3 


HIP062145 


K3V 


1.50 


94 


1.31 


4.84 


90 


0.91 


0.96 


1.31 


4.66 




10.12^ 


0.34 


15.53 


N 


2 


1;2;6 


HIP062207 


GOV 


6.40 


117 


13.68 


6.21 


116 


1.17 


0.99 


58.50 


2.49 


18.00 


47.26 


3.16 


12.99 


Y 


3 


1;12;15;23 


HIP062523 


G5V 


1.00 


101 


11.14 


5.37 


99 


1.00 


0.98 


20.40 


2.35 


3.94 


3.57 


2.72 


-0.66 


N 


3 


2;9;10;11;12;15 


HIP064792 


GOV 


0.34 


219 


25.64 


11.47 


226 


2.27 


1.03 


13.60 


5.06 


-2.38 


11.94 


3.92 


0.12 


N 


3 


1;2;4;5;9;12;13;23 


HIP064797 


K1V+M1V 


1.00 


128 


14.52 


6.66 


128 


1.29 


1.00 


9.61 


5.40 


-0.91 


5.32^ 


3.12 


-0.43 


N 


3 


1;2;10;13;15;23 


HIP065026 


KO 


0.80 


143 


16.34 


7.85 


149 


3.76 


1.04 


- 






14.56^ 


3.75 


1.79 


N 




2 


HIP065721 


G5V 


8.30 


382 


43.72 


19.74 


386 


3.86 


1.01 


41.20 


3.15 


-0.80 


35.11 


2.72 


5.65 


Y 


3 


1;4;6;9;10;12;23;24 


HIP067275 


F6IV+M2 


1.30 


340 


38.52 


18.21 


340 


3.40 


1.00 


32.70 


5.54 


-1.05 


0.86 


5.34 


-3.25 


N 


3 


1;2;4;9;10;13;16 


HIP067422 


K4V+K6V 


0.85 


101 


12.93 


6.27 


101 


1.03 


1.00 


- 






6.77^ 


2.79 


0.18 


N 


3 


1;2;15 


HIP067620 


G5+V 


2.30 


96 


11.42 


5.38 


101 


0.70 


1.05 


8.97 


2.75 


-0.89 


4.77^ 


3.40 


-0.18 


N 


3 


2;3;5;9;10 


HIP068184 


K3V 


5.50 


170 


19.61 


9.09 


173 


1.73 


1.02 


17.70 


6.59 


-0.29 


8.71 


2.96 


-0.13 


N 


2 


2;9 


HIP068682 


G8V 


5.00 


123 


14.01 


6.62 


124 


1.24 


1.01 


4.62 


5.49 


-1.71 


6.93^ 


3.06 


0.10 


N 


3 


1;5;6;12;15 


HIP069965 


F9V 


1.50 


115 


13.15 


6.44 


116 


1.28 


1.00 


- 






10.46^ 


3.22 


1.25 


N 


1 


3 


HIP070319 


G1V 


5.20 


99 


11.06 


5.23 


98 


0.99 


0.99 


5.21 


6.03 


-0.97 


5.20^ 


2.79 


-0.01 


N 


3 


1;5;9;12;15;24 


HIP070857 


G5 


3.60 


79 


-16.16 


4.20 


79 


0.79 


0.99 


6.87 


4.49 


5.13 


5.16^ 


2.83 


0.34 


N 


2 


9;12 


HIP071181 


K3V 


3.00 


93 


9.83 


4.68 


86 


0.86 


0.92 


33.80 


4.70 


5.10 


13.94^ 


2.51 


3.69 


Y 


2 


1;6 


HIP071681 


K1V 


6.00 


39085 


4472.81 


2008.32 






- 


- 




-0.45 


1809.12 


33.20 


-6.00 


N 


3 


2;3;4;5;23;28 


HIP071908 


A7V 




608 


71.05 


34.50 


632 


6.32 


1.04 


- 






36.94 


5.96 


0.41 


N 






HIP072567 


G1V 


0.45 


119 


13.42 


6.34 


119 


1.20 


1.00 


11.10 


2.94 


-0.79 


5.44 


2.80 


-0.32 


N 


3 


2;5;9;10;11;12 


HIP072603 


F3V 


0.30 


142 


17.72 


8.69 


143 


1.44 


1.01 


- 






5.58^ 


2.59 


-1.20 


N 


3 


2;18;29 


HIP072848 


K2V 


0.40 


199 


22.94 


10.66 


201 


2.01 


1.01 


35.80 


5.59 


2.30 


20.02 


2.69 


3.48 


Y 


3 


2;12;15;23 


HIP073100 


F7V 


4.00 


126 


14.83 


6.68 


131 


1.32 


1.04 


25.80 


2.11 


5.20 


14.39 


3.12 


2.47 


Y 


3 


2;4;15;24 


HIP073184 


K4V 


1.10 


452 


47.61 


2.17 


420 


4.20 


0.93 


52.70 


3.18 


1.60 


16.30 


4.18 


3.38 


N 


3 


2;3;5 


HIP073996 


F5V 


1.80 


208 


23.59 


10.81 


208 


2.08 


1.00 


36.70 


4.40 


2.98 


13.60 


4.23 


0.66 


N 


3 


1;2;16;23 


HIP077052 


G2.5V 


3.00 


148 


16.97 


8.39 


145 


5.31 


0.98 


- 






2.06 


3.33 


-1.90 


N 


3 


2;4;5;9;10;28 


HIP078459 


GOV 


7.70 


204 


23.10 


10.67 


204 


2.04 


1.00 


30.70 


4.13 


1.84 


1.86 


4.76 


-1.85 


N 


3 


1;9;10;12;15;23;24 


HIP078775 


G8V 


6.50 


92 


10.12 


4.80 


89 


0.89 


0.97 


13.00 


2.85 


1.01 


-8.82^ 


3.71 


-3.67 


N 


3 


1;6;9;10;12;15;23 


HIP079248 


KOV 


8.20 


96 


10.63 


5.05 


94 


0.95 


0.98 


10.30 


1.85 


-0.18 


7.36 


2.85 


0.81 


N 


3 


1;6;9;10;15 


HIP080725 


K2V 


1.50 


77 


8.80 


4.40 


74 


3.89 


0.96 


- 






4.56^ 


3.29 


0.05 


N 


1 


2 


HIP082860 


F8V 


1.70 


276 


30.93 


14.63 


273 


2.73 


0.99 


47.20 


4.93 


3.30 


43.51 


2.75 


10.50 


Y 


3 


1;2;23 


HIP083389 


G8V 


4.30 


72 


7.97 


3.83 


71 


0.71 


0.98 


5.33 


5.08 


-0.52 


7.29^ 


2.66 


1.30 


N 


3 


1;6;9;12 


HIP084862 


GOV 


6.90 


220 


24.62 


11.58 


218 


2.18 


0.99 


23.90 


4.25 


-0.17 


8.95 


4.46 


-0.59 


N 


3 


1;9;12;24 


HIP085235 


KOV 


5.60 


116 


13.55 


6.08 


113 


1.14 


0.97 


53.80 


1.61 


25.00 


29.95 


2.75 


8.68 


Y 


2 


5;12 


HIP085295 


K7V 


1.10 


176 


19.48 


9.17 


172 


1.73 


0.98 


18.30 


3.02 


-0.39 


1.58 


5.75 


-1.32 


N 


3 


1;2;15 


HIP086036 


GOVa 


1.00 


253 


30.40 


13.50 


253 


2.53 


1.00 


30.40 


4.62 




13.54 


4.45 


0.01 


N 




2 


HIP086796 


G3IV-V 


7.70 


267 


30.77 


15.96 


272 


2.72 


1.02 


31.50 


7.29 


0.10 


13.59 


5.38 


-0.44 


N 


3 


3;5;24 


HIP088601 t 


KOV 


1.20 


879 


106.34 


79.78 


870 


8.70 


0.99 


126.00 


5.40 


3.64 


55.46 


5.75 


-4.23 


N 


3 


1;2;4;5;15;23 


HIP088972 


K2V 


5.80 


149 


16.95 


7.95 


150 


1.51 


1.01 


9.30 


4.58 


-1.67 


7.00^ 


3.28 


-0.29 


N 


3 


1;2;5;6;10;13;15 


HIP089042 


GOV 


5.10 


174 


19.86 


10.18 


160 


8.84 


0.92 








9.75 


3.56 


-0.12 


N 


3 


3;5;24 


HIP091009 


K6Ve 


0.01 


88 


9.93 


4.87 


88 


0.88 


1.00 


7.68 


2.42 


-0.93 


4.62^ 


2.77 


-0.09 


N 


3 


2;15 


HIP092043 


F6V 


2.20 


455 


49.93 


24.57 


441 


4.41 


0.97 


69.30 


8.46 


2.29 


28.35^ 


8.22 


0.46 


N 


3 


1;2;4;5;16;23 


HIP095995 


K2V 


7.50 


122 


12.82 


6.31 


113 


1.13 


0.93 


13.60 


5.20 


0.15 


1.29^ 


3.04 


-1.65 


N 


1 


15 


HIP096100 


G9V 


3.50 


591 


66.96 


29.81 


591 


5.91 


1.00 


74.70 


5.82 


1.33 


28.80 


5.94 


-0.17 


N 


3 


1;2;9;10;12;15;23 
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Table 2 — Continued 



o 



MIPS 



PACS 



Far 



Name 


Sp. 


Age 


P 2 4 


P70 


P100 


F24 


°24* 


R24 


^70 


070 


X70 


^100 


^100 


X100 


IR 


Age 


Age ref. 




type 


(Gyr) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


Exc? 


Flag 




t nnnn/c a a 1 

H1P096441 


F4V 


1.50 


290 


33.86 


15.91 


299 


3.00 


1.03 


40.30 


2.73 


2.36 


15.43 


5.35 


-0.09 


N 


1 


2;23;24 


111 T}f\(\^l(\ A A 

H1P097944 


K3V 


0.75 


197 


21.75 


10.70 


201 


2.01 


1.02 








12.14 


4.36 


0.33 


N 


1 


2;3 


HIP098959 


G2V 


4.50 


116 


12.94 


6.16 


115 


1.16 


0.99 


23.50 


2.64 


4.00 


13.68^ 


3.08 


2.44 


Y 


3 


3;5;24 


HIP099461 


K2.5V 


7.00 


481 


52.92 


29.11 


467 


4.67 


0.97 


50.70 


8.52 


-0.26 


24.62 


5.75 


-0.78 


N 


3 


2;3;5;12;17;23 


H1P101955 


K5V 


0.90 


106 


11.30 


5.47 


102 


1.03 


0.96 


16.60 


6.02 


0.88 








N 


1 


2 


H1P101997 


G8V 


5.50 


106 


11.68 


5.52 


103 


1.05 


0.97 


3.73 


2.83 


-2.81 


2.76 


2.85 


-0.97 


N 


3 


3;4;5;9;12 


HIP 103 3 89 


F6V 


0.45 


118 


13.80 


6.24 


122 


1.23 


1.03 


45.30 


2.10 


15.00 


23. 17^ 


2.85 


5.94 


Y 


3 


2;3;12;16 


HIP104214^ 


K5V 


1.10 


932 


96.25 


45.34 


885 


8.85 


0.95 








37.54 


6.29 


-1.24 


N 


2 


1;2;5;10;15 


HIP105312 


G7V 


6.40 


112 


12.81 


5.94 


113 


1.14 


1.01 


11.00 


6.72 


-0.27 


9.72^ 


3.10 


1.22 


N 


1 


3 


HIP 106696 


K1V 


1.80 


77 


8.22 


3.88 


73 


0.74 


0.95 


11.50 


5.29 


0.62 


7.14^ 


3.33 


0.98 


N 


3 


2;3;5;17 


HIP107350 


GOV 


0.35 


111 


12.96 


5.86 


114 


1.15 


1.03 


25.50 


2.56 


4.90 


9.79 


2.71 


1.45 


Y 


3 


1;2;3;5;12;13;15;16;23 


HIP107649 


G2V 


4.00 


161 


19.46 


8.50 


166 


1.67 


1.03 


398.00 


7.01 


54.00 


236.22^ 


3.58 


63.61 


Y 


3 


2;3;5;17 


HIP108870 


K5V 


2.00 


1147 


123.29 


53.88 


1090 


10.90 


0.95 


112.00 


6.07 


-1.86 


53.96 


2.71 


0.03 


N 


3 


2;5;15;17 


HIP109378 


GO 


8.10 


87 


9.64 


4.48 


85 


0.86 


0.97 


10.30 


2.00 


0.33 


7.24 


2.82 


0.98 


N 


3 


1;9;10;12;15 


HIP109422 


F6V 


4.90 


220 


24.94 


15.64 


220 


2.21 


1.00 


8.61 


5.69 


-2.87 


13.50 


5.36 


-0.40 


N 


3 


3;9;16;24 


HIP1 13576 


K7+Vk 


1.10 


139 


21.70 


6.97 


135 


1.36 


0.97 


21.70 


2.38 




5.53 


3.34 


-0.43 


N 


1 


2 


HIP1 14948 


F6V 


0.33 


119 


10.64 


6.69 


141 


1.42 


1.18 


72.70 


1.83 


33.91 


36.89 


2.86 


10.56 


Y 


2 


2;3 


HIP1 16745 


K3+V 


3.50 


113 


12.11 


5.57 


107 


1.08 


0.95 


16.20 


1.77 


2.31 


4.59 


2.87 


-0.34 


N 


2 


2;3;17 


HIP 120005 


M0.0V 


0.44 


374 


42.75 


22.97 














17.04 


3.47 


-1.71 


N 


1 


2 



References. — (1) Duncan et al. 1 1991); (2) Rosat All Sky Survey; (3) Grav et al. 1 2006); (4) Schroder et al. 1 2009); (5) Henrv et al. 1 1996); (6) Rocha-Pinto & Made] 1 1998); (7) Vican 1 2012) - isochrone ages; (8) Schmitt & Liefke O 
2004); (9) Wright et al. 1 2004); (10) Katsova & Livshits 1 201 1); (11) Martmez-Arnaiz et al. 1 2010); (12) Isaacson & Fischer 1 2010); (13) Vican 1 2012) - gyro ages; (14) Barnes 1 2007); (15) Grav et al. 1 2003); (16) v sin(z); (17) Jenkins et al, 



120061); (T8)lMontes et alJ feOOlh; Q9)IVicanU2012j) - X- rav; (20) White et al. 12007); (21) log(g); (22) Lachaume et al. 1 1999); (23) Buccino & Mauas 12008); (24) HR diagram position; (25) /3 Pic MG; (26) Nakaiima et al. i 201(3); (27) 
Jenkins et al. < 2011); (28) Mamaiek & Hillenbrand i 2008); (29) Barrado y Navascues 1 1998) 

* HD007439: Used standard colors in place of K.; HD061421: K magnitude derived from COBE measurements. 

^ K band data used instead of W3. 



PACS 70 /im data also available and was used for Spitzer MIPS 70 /im comparison. 

HD109085: PACS 70 fim data: F 70 = 186.89 mJy, a 70 = 6.95 mJy, X io =21.16 mJy; HIP007978: PACS 70 /im data: F 1Q 
When nominal error in the far-IR was < 5%, we assumed a 5% error. 
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= 735.09 mJy, a 10 = 6.53 mJy, xio = 110.06 mJy 



The age flag gives the number of independent methods yielding consistent age values. 
** Allowing for systematics, 1% photometric error was root-sum-squared with the statistic ones at 24 /xm lEngelbracht et al. 2007). 



